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ABSTRACT 

We present comprehensive photometric and spectroscopic observations of the faint transient 
SN 2008S discovered in the nearby galaxy NGC 6946. SN 2008S exhibited slow photometric 
evolution and almost no spectral variability during the first nine months, implying a long 
photon diffusion time and a high density circumstellar medium. Its bolometric luminosity ( 
~ 10 41 erg s _1 at peak) is low with respect to most core collapse supernovae but is comparable 
to the faintest type II-P events. Our quasi-bolometric lightcurve extends to 300 days and shows 
a tail phase decay rate consistent with that of 56 Co. We propose that this is evidence for an 
explosion and formation of 56 Ni (0.0014 ± 0.0003M Q ). Spectra of SN 2008S show intense 
emission lines of Ha, [Ca II] doublet and Ca II NIR triplet, all without obvious P-Cygni 
absorption troughs. The large mid-infrared (MIR) flux detected shortly after explosion can be 
explained by a light echo from pre-exisiting dust. The late near-infrared (NIR) flux excess is 
plausibly due to a combination of warm newly-formed ejecta dust together with shock-heated 
dust in the circumstellar environment. We reassess the progenitor object detected previously 
in Spitzer archive images, supplementing this discussion with a model of the MIR spectral 
energy distribution. This supports the idea of a dusty, optically thick shell around SN 2008S 
with an inner radius of nearly 90AU and outer radius of 450AU, and an inferred heating source 
of 3000 K. The luminosity of the central star is L ~ 10 46 L©. All the nearby progenitor dust 
was likely evaporated in the explosion leaving only the much older dust lying further out in the 
circumstellar environment. The combination of our long term multi-wavelength monitoring 
data and the evidence from the progenitor analysis leads us to support the scenario of a weak 
electron capture supernova explosion in a super-asymptotic giant branch (AGB) progenitor 
star (of initial mass 6-8 M©) embedded within a thick circumstellar gaseous envelope. We 
suggest that all of main properties of the electron capture SN phenomenon are observed in SN 
2008S and future observations may allow a definitive answer. 

Key words: supernovae: general - supernovae:individual: SN 2008S - super- 
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1 INTRODUCTION 



Table 1. Properties of NGC6946. 



In recent years deeper and more frequent searches for transient 
events and stellar explosions in the local and distant Universe have 
provided us with important information on the evolution of the 
most massive stars. However, the simplicity of the emerging pic- 
ture is compromised by the growing number of peculiar even ts (e.g. 
iKulkarni et alfcOOllSmith et alfeOOllOuimbv et alj|2007l) . 

From an observational point of view, the challenge is to de- 
cide when the introduction of new classes is required or if peculiar 
or novel transients are just variations of an understood scheme. The 
discovery of some low energy events (in terms of their bolometric 
luminosity and kinetic energies) leads us to investigate in more de- 
tail the observational differences between explosive (core collapse, 
pair instability explosions) and eruptive (pair instability pulsation, 
outburst) transients. 

From a theoretical point of view, recent observations demon- 
strate that the standard scenario of stellar evolution and explosion 
physics may not be complete. Both the ex tremely bright type II SNe 
(Woos lev et"ai]|2q07l:lLanger et alj|2007t) and the faint t ype II SNe 
( Smith et al Jl2008l ; lBond et alJl2009l ; lBerger et alj|2009h have been 
proposed to have physical origins other than the core-collapse of a 
degenerate Fe (or O-Ne-Mg) core. 

SN 2008S is one of the most intriguing transient events dis- 
covered in recent years. Although it has been given a supernova 
designation, (which we will employ in this paper) it is not yet 
certain that it was a supernova of the canonical core-collapse 
type (CCSN). The tra nsient was discovered in NGC 6946 by 
I Arbour & Bolesl d2008l) on February 1.78 UT with a 30-cm f/6.3 
Schmidt-Cassegrain reflector at about 17.6 mag. Eight confirming 
images of SN 2008S were taken on February 2 .76 UT, yielding 
a magnitude of 17.1. Furthermore, I Arbour! b008h provided a new 
ima ge of SN 2008S acquired on 2008 January 24 UT (17.8 mag) 
and lSchmeed j2008h reported an image obtained on 2008 January 
30.529 (16.7 mag). T he transient wa s class ified as a young red- 
dened Type Iln SN by Sta nishev et al] J2008I) based on a low reso- 
lution spectrum taken at the Nordic Optical Telescope, with narrow 
H/3 an d Ha emission lines and strong Na I D doublet. Steele et al. 
( 2008) reported a new spectrum of SN 2008S obtained on Feb. 29 
UT with the 3-m Shane reflector equipped with Kast double spec- 
trograph at the Lick Observatory, and suggested SN 2008S to be a 
"SN impostor" based on peculiar spectral properties and the very 
faint absolute visual magnitude. 

Remarkably, a bright point-like source coincid ent with SN 
2008S was detected in archival Spitzer MIR images bv lPrieto et al.l 
d2008l) . They found no optical counterpart to this precursor and 
suggested this MIR source was a stellar progenitor with mass of 
about 10 M© and luminosity of ~ 3.5 x 1O 4 L0, enshrouded 
in its own dust. The stellar mass and the total luminosity esti- 
mates result from a blackbody fit to the MIR spectral energy dis- 
tribution (SED) of the progenitor star. Shortly afterwards, another 
transient was discovered in the nearby galax y NGC 300 which 
bears a striking r esemblance to SN 200 8S teerger et alj 1 20091 ; 
Bond et al. 2009]). IThompson etall (2008) reported the discovery 
of a similar progenitor star in Spitzer MIR prediscovery images 
and again an optical counterpart was lacking. They suggested that 
both transients share a common evolutionary channel and also 
that t he optical transient discovered in M 85 was of similar ori- 
in jKulkami et alj|2007l; IPastorello et alj|2007l) . IThompson et al.l 
2008) and IPrieto et al] ( l2008l) have proposed that these events 
could be low energy electron-capture SNe (ECSNe) from stars of 
initial mass around 9M©. The existence of such explosions has 
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a based on the local velocity field model given in Mould et al. (2000) using 
the terms for the influence of the Virgo Cluster, the Great Attractor, and the 
Shapley Supercluster. 



been theoretica lly predicted for many years dMivaii et al 



Nomoto 1984; Mivaii & Nomotol 119871; lHashimoto et al 
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iKitaura et a l. 2006; Po elarends et al.l20 08). However, what exactly 
the mass range of the progenitors would be and how the SN evo- 
lution would appear is far from certain. The nature of these tran- 
sients has not yet been firm ly established, since rece nt works on 
SN 2008S dSmith et alj2008h and NGC 300 OT2008-1 dBond et al.l 
2009; Berger et al. 2009) suggest that these events are the outbursts 
of a massive star and not the cataclysmic stellar deaths of stars after 
core-collapse. 

In this paper we present results from our extensive photomet- 
ric and spectroscopic follow-up of SN 2008S, together with anal- 
ysis of supernova and progenitor observations. The properties of 
the host galaxy are described in Sect. [2] Photometric data reduction 
and analysis are detailed in Sect. [3] and the evolution of the SED 
is illustrated in Sect. [4] Spectroscopic data reduction and analysis 
are detailed in Sect. [5] Sect.|6]is devoted to the analysis of the pre- 
explosion images and to the discussion of the progenitor star. A 
summary of our observations, a comparison with other underlumi- 
nous transients and some type II-L SNe, and our conclusions on the 
nature of SN 2008S are given in Sect.|7] 



2 HOST GALAXY, DISTANCE AND EXTINCTION 

SN 2008S was discovered at R.A.= 20 h 34 m 45 i :37 and Dec.= 
60°05'58"3 (2000), about 53 arcsec West and 196 arcsec South 
of the nucleus of NGC 6946. Details of the host galaxy obtained 
from the NASA/IPAC Extragalactic databasfl are summarised in 
TableCO 

Optical, far infrared, radio continuum and X-ray observations 
indicate vigorous star formation (SF) throughout the NGC 6946 
disc, one of the highest among nearby spiral galaxies, a mild star- 
burst at its centre, and an interstellar medium (ISM) stirred by SNe 
and stellar winds l|En gargiola 19 9ll;|Boulanger & Viallefonol 19921 ; 
iKamphuis & Sancisill993l ; ISchlegelll994l ; lLacev et al.lll997l) . This 
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© 0000 RAS, MNRAS 000, 000-000 



SN2008S 



high level of SF in the disc of NG C 6946 has been attrib uted 
both to its strong spiral density wave (Tacconi & Young 1990) and 
to stochastic, self-propagating SF dDegioia-Eastwood et al.ll 19841) . 
Signs of low level SF such as HII regions and UV bright clusters 
have been discovered in the far outer regions of galactic disc well 
beyond the R25 radius. 

Eight other SNe have been detected in this galaxy six of which 
were classified as type II SNe (1917A (II), 1948B (II-P), 1968D (II) 
1980K (II-L) 2002hh (II-P) 2004et (II-P)) and two remain unclas- 
sified (SN 1939C and SN 1969P). All these SNe were brighter than 
mag 15 except SN 2002hh which was highly reddened. Among 
these eight SNe, four have been detected as radio SNe (1968D, 
1980K, 2002hh, 2004et) and three as X-ray SNe (1968D, 1980K, 
2004et). Many SN remnants have been detected in NGC 6946 us- 
ing optical, radio and X-ray telescopes dMatonick & Feseijl 19971 ; 
ISchlegel et alj|200d : |pannuti et al.ll2007h . 



2.1 Metallicity 

The galactic metallicity at the position of SN 2008S can be esti- 
mated in a similar way to that for two other recent SNe in this 
galaxy (SN 2002hh and SN 2004et) as shown by Sma rtt et al 



2009). The abundance gradient determined by Pilvug in et al 
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20041) (12 + logQ/H = 8.7-0.41(i?/i? 25 )) and the de-projected 
galactocentric radius of the SN position can be used to determine 
the likely local metallicity at the position of SN 2008S. Using the 
distance of 5.7 Mpc as discussed in Sect. 12.21 SN 2008S is at a de- 
projected galactocentr ic radius of 4.9 kp c, and with R25 = 9.1 kpc 
(fr om HyperLedfl IPaturel et alj d2003l) ). the metallicity gradient 
of IPilvugin et alj ( 120041) results in an appro ximate oxygen abun- 
dance of 8.5 dex. On the lPilvugin et alj ( 2004) abundance scale, so- 
lar is approximately 8.7 dex. Hence the environment of SN 2008S 
is mildly sub-solar, although within the uncertainties in this method 
a solar-like composition of the progenitor is still quite possible. 
By comparison, the oxygen abundances estimated for SN 2002hh 
and SN 2004et ar e approximately 8.5 and 8.3 dex, respectively 
dSmartt et alj|2009h . 



2.2 Distance 

There are several estimates of the distance to NGC 6946, obtained 
with different methods and listed in Table [2] Two S Ne hosted in 
NGC 6946 have been used as distance indicators. ISchmidt et alj 
(1994) applied the Expanding Photosphere Method (EPM) to 
SN 1980K and f ound a distance modulus of 28.78 ± 0.4 mag. 
ISahuet alj (2006) used a "standard candle method" (SCM) for SN 
2004et, based on the correlation between the expansion velocities 
of the SN II- P ejecta and the bolometric luminosities during the 
plateau phase dHamuv & Pintol2002l : lNugent et al.l2006h . They ob- 
tained fi — 28.78 ± 0.11 mag, in close agreement with the esti- 
mate obtained with EPM for SN 1980K. There is one distance esti- 
mate that is significantly different from the rest: the radio observa- 
tions of SN 1980K yield a much larger value of 30.5 ± 0.3 Mpc 
( Weile ret alj 1 1998h . As this is much larger than the other esti- 
mates (which are consistent within the uncertainties) we shall dis- 
count this value and use an unweighted mean of 28.78 ± 0.08 mag 
throughout this paper. 



20 40 60 
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Figure 1. Temporal evolution of the EW(Na I D). Phase is in days after the 
explosion epoch (JD 2 454 486). 



Table 3. Measurements of the EW of Na I D. 



JD a 


ph 6 


EW (A) 


501 


16 


4.4 ±0.50 


504 


19 


3.4 ±0.45 


508 


22 


3.2 ±0.45 


514 


28 


2.7 ±0.35 


524 


38 


2.5 ±0.48 


530 


44 


2.5 ±0.40 


533 


47 


2.3 ±0.48 


554 


68 


2.6 ±0.48 



2 HyperLeda database, http://leda.univ-lyonl.fr 
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a JD- 2454000.00 

b Phase is in days after the explosion date JD 2 454 486 ± 4. 



2.3 Extinction 

NGC 6946 is located close to the Galactic plane (Galactic latitude 
~ 12°), with an estim ated reddening of E(B — V) = 0.342 mag 
dSchlegel et al.|[l99c3) . For the SNe which occurred in this galaxy, 
different values of extinction have been estimated depending on the 
SN position. In all cases the presence of the Na I D (AA5890,5896) 
lines has been used as an indicator of the prese nce of dust and used 
to estimate the reddening at the SN position 1 Zwitter et alj|2004l 
ISahu et al.ll200o1 : lMeikle et alJl2006l:|Pozzo et alj|200rj) . Strong Na 
I D lines in absorption are also present in the SN 2008S spectra 
until about 70 days after the explosion. The EW(Na I D) appears 
to show a temporal evolution from 4.4 A to 2.5 A during this time 
(Tableland Fig.Q}. 

To estimate the error in our EW measurements we performed 
a Monte-Carlo simulation adding a number of absorption lines with 
known EW at different positions in each SN spectrum and re- 
measured their EW. We repeated the simulation for different val- 
ues of line EW and strength. These simulations were performed 
separately in each spectrum to take account of the differences in 
spectral resolution and S/N ratio. In order to investigate the reality 
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Table 2. Estimates of the distance to NGC 6946. 



rv 

Distance 


TV A 

Distance modulus 


ivieinou 


_ 

Reference 


(mpcj 


(mag) 






5.5 


28.70 


HI Tully-Fisher relation 


Pierce (1994) 


5.4 


28.66 


CO Tully-Fisher relation 


Schonieer & Sofue (1994) 


6.0 ±0.5 


28.90 ±0.18 


galaxy brightest supergiants 


Sharinaetal. (1997) 


5.9 ± 0.4 


28.85 ±0.15 


brightest supergiants of group 


Karachentsev et al. (2000) 


5.6 ± 1.8 


28.73 ±0.68 


"sosie galaxies" 


Terry et al. (2002) 


6.1 ± 0.6 


28.92 ±0.21 


planetary nebulae luminosity function 


Herrmann et al. (2008) 


5.7 ±0.7 


28.78 ±0.40 


EMP SN 1980K 


Schmidt et al. (1994) 


5.7 ±0.3 


28.78 ±0.11 


SCM SN 2004et 


Sahu et al. (2006) 



of the EW changes in the Na I D lines, we carried out a quantitative 
statistical test, performing a linear fit to the data, and found a neg- 
ative slope at the 2 sigma (95% confidence) level. To test further 
that the data are better represented by a temporally declining EW, 
rather than a fixed value, we exploited the Bayesian information 
criteri on (BIC), w hich give an approximation for the Bayes factor 
(see lLiddldl2004 and references therein). The BIC is defined as 
X 2 + Npar log Ndata where x 2 is the total x 2 for the model, N par 
is the number of parameters of the model and Ndata is the number 
of data points used in the fit. The best model minimizes the BIC. A 
difference of 2 for the BIC is regarded as positive evidence, and of 6 
or more as strong evidence, against the model with the larger value. 
The BIC corresponding to no evolution is larger by 6 than the BIC 
for a straight line with slope (ABIC = BIC CO nst — BIC s i opc > 6), 
clearly supporting the scenario for a decreasing temporal evolution. 
The two BIC values are comparable only if we exclude from the 
analysis the first two epochs. We conclude that we find evidence 
for a change in EW of the Na I D feature. 

The evolution of the local component EW(Na I D) may be due 
to an evolution of the ionization conditions in the CSM and in the 
ejecta of SN 2008S since the EW is related to the ionization stage 
of Na I. The evolution in the EW may also imply that the local 
extinction underwent a temporal decline. However, given (a) the 
lack of any well-established EW(Na I D)-exinction correlation at 
the very large EWs involved, and (b) the possibility that the Na I D 
feature includes saturated components, we make no attempt to use 
the EW(Na I D) to determine the extinction or its possible variation. 
In any case, the EW(Na I D) variation could be simply due to evo- 
lution of the physical properties of the gas around SN 2008S (see 
Sect. |5.2| l, with the extinction taking place at a completely different 
location. By days 182 and 256 the Na I D feature has become visi- 
ble in emission. This change is also indicative of the circumstellar 
origin of this feature. 

The presence of circumstellar Na I D has also been observed 
in the type Iln SN 1998S and has been interpreted as a signature 
of slow moving outflows originating from the progenitor while its 
blueshift and growing intensity between 20-40 days after explo- 
sion h as been associated w ith variable physic a l cond itions in the 
CSM ( lBowenetal.1l2000l) . IChugai & Utrobinl J2008l) studied the 
formation of the Na I D and Ca II H&K lines in the RSG wind 
after a SN II-P explosion with the goal of using these as a diag- 
nostic of the wind density. They extrapolated their model to a very 
high wind density to reproduce the intensity of these lines observed 
in SN 1998S. However, the EW of the absorption depends non- 
monotonically on the wind density. The case of SN 1998S with its 
very dense wind has shown that the EW(Na I D) decreases with 
wind density because of the ionization of metals in the wind by 
UV radiation. Variable EW(NaI D) have been detected also in a 



few Ty pe la SNe (2006X. |patat et al.1 d2007t) : 1999cl. lBlondin et all 
(2008); 20071e, Simon et al. in pre p). However, this do es not seem 
to be a very common phenomenon dBlondin et alj 2008) and the in- 
terpretation in terms of evolution of the CSM physi cal conditions 
indu ced by the SN radiation field jPatat et alj|2007l) is still debat- 
able dChugaill2008h . 

Here we adopt the Galactic absorption in V band, Ay = 1.13 
mag, calculated from the list of A/E(B — V) of Schle geTet al.l 
( 1998) along with their estimate of E(B — V) and an extinction 
local to the SN, Av ~ 1 mag, required by our light echo model to fit 
the observed SED at 17 days after explosion when the MIR excess 
is observed, (see Sect. 14. U . If we assume Av = 2.13 mag, the EW 
of K I (A 7699) wou l d be a bout 0.17 A following the calibration 
bv lMunari & Zwitten 1 d 19971) . Unfortunately, the K I region of the 
spectrum lies close to strong telluric absorption (7570-7750 A) and 
is only ever covered at low resolution so we did not observe the K 
I feature. 



3 PHOTOMETRIC DATA AND ANALYSIS 

We commenced monitoring SN 2008S shortly after the discovery 
epoch and collected data for the following eight months with a sam- 
pling rate among the highest ever obtained for such a peculiar tran- 
sient. Data obtained before the discovery date by several amateur 
astronomers are also included to constrain the explosion epoch of 
SN 2008S. The unfiltered image acquired by D. Abraham on Jan- 
uary 16 (JD 2 454 482) shows no object visible in the SN 2008S 
location with a limiting magnitude of 19.20 in V band (18.60 in 
R band), while the first detection of SN 2008S is eight days later 
on January 24 (JD 2 454 490). We therefore adopt January 20 (JD 

2 454 486) as the explosion epoch, the uncertainty being about 4 d. 
The phases in this paper are relative to the explosion date (when we 
fix ph=0). 

3.1 Optical data 

Optical photometry of SN 2008S was obtained with many tele- 
scopes and a summary of their characteristics is given in Table [4] 
Unfiltered images were obtained for many epochs with a 40 cm 
telescope with a SXVF H9 camera, a 30 cm telescope with a 
MX916 camera, a 35 cm telescope with a ST-9E/9XE camera and 
with a 25 cm telescope with a ICX424 CCD. 

Basic data reduction (overscan correction, bias subtraction, 
flat fielding, trimming) was performed using standard routines 
in IRAlfl The instrumental magnitudes were obtained with the 

3 Image Reduction and Analysis Facility (IRAF) is distributed by the Na- 
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Table 4. Summary of the characteristics of the telescopes used during the photometric follow up. 



Telescope 


Primary mirror 


Camera 


array 


CCD 


pixel scale 


field of view 


filters 




m 








arcsec/pix 


arcmin 




TNGD 


3.6 


DOLORES 


2048 x 2048 


EEV 42-40 


0.25 


8.6 x 8.6 


Johnson U, B, V; Cousin R, I 


TNGN 


3.6 


NICS 


1024 X 1024 


HgCdTe Hawaii 


0.25 


4.2 x 4.2 


JHK 


NOT 


2.5 


ALFOSC 


2048 x 2048 


EEV 42-40 


0.19 


6.4 x 6.4 


Johnson U,B,V,R 


CAHAT 


2.2 


CAFOS 


2048 x 2048 


SITe 


0.53 


16 x 16 


Johnson B, V, R, I 


LT 


2.0 


RATCam 


2048 x 2048 


EEV 42-40 


0.13 


4.6 


Bessel U, B, V; Sloan r',i' 


CAO 


1.8 


AFOSC 


1024 X 1024 


TK1024AB 


0.46 


7.8 


Bessel B,V,R\ Gunn i 


SAO 


1.0 




2048 x 2048 


EEV 42-40 


0.48 


8.3 


Johnson V 


AZT24 


1.0 


SWIRCAM 


256 x 256 


HgCdTe PICNIC 


1.03 


4.4 


JHK 


MSK 


0.7 


Apogee AP-7p 


512 x 512 


SITe 


0.94 


4 


Johnson B, V, I; Cousin R 


MSKL 


0.7 


Apogee AP-7p 


512 x 512 


SITe 


0.64 


5.5 


Johnson B, V, I; Cousin R 


CRM 


0.6 


Apogee AP-47p 


1024 x 1024 


Marconi47-10 


0.71 


6.1 


Johnson B, V, I; Cousin R 


SLV 


0.5 


SBIG ST-10XME 


2184 X 1472 


KAF3200ME 


1.12 


20.6 x 13.9 


Johnson B, V; Cousin R, I 


WOT 


0.4 


SBIG ST-10XME 


2184 x 1472 


KAF3200ME 


0.44 


16.0x10.8 


Sloan r' 



TNGD = the Telescopio Nazionale Galileo (TNG) with the Device Optimized for the LOw RESolution (DOLORES); TNGN = the Telescopio Nazionale 

Galileo (TNG) with the Near Infrared Camera Spectrometer (NICS); NOT= the Nordic Optical Telescope (NOT) with the Andalucia Faint Object 
Spectrograph and Camera (ALFOSC); CAHAT= the 2.2 m telescope at Calar Alto Observatory (CAHA) with the Calar Alto Faint Object Spectrograph 
(CAFOS); LT = the Liverpool Telescope (LT) with the optical CCD camera RATCam; CAO = the Copernico telescope at Asiago Observatory with the 
Asiago Faint Object Spectrograph and Camera (AFOSC); SAO = the 1 m telescope of Special Astrophysical Observatory of Russian Academy of Sciences; 
AZT24= the AZT 24 telescope at Campo Imperatore Observatory with SWIRCAM; MSK = the 70 cm telescope of the Sternberg Astronomical Institute in 

Moscow; MSKL = the 70 cm telescope of the Sternberg Astronomical Institute +focal reducing lens; CRM = the 60 cm reflector of the Sternberg 
Astronomical Institute Crimean laboratory; SLV = the 50 cm telescope of the Astronomical Institute of Slovak Academy of Sciences at Tatranska Lomnica; 

WOT = the 40 cm telescope at the Wendelstein Observatory 

point spread function (PSF) fitting technique using a custom made 
DAOPHOT based package (SNOoPY). We did not apply the tem- 
plate subtraction technique, since the host galaxy contamination is 
negligible around SN 2008S in the optical and in the NIR range. 

The photometric calibration was carried out by a comparison 
with Landolt standard stars observed the same night when possible. 
As our local sequ ence of stars, we chose a subset of that adopted in 
IPozzo et al.l feOOfj) for SN 2002hh, shown in Fig.[2]and calibrated 
it with respect to a number of Landolt standard fields on several 
photometric nights. 

Our calibration is in agreement with that o f lPozzo et al.N200fj|) 
with an average difference of 0.01 mag in the V and I bands and of 
0.005 mag in the 7? band. The magnitudes of the calibrated local se- 
quence, listed in Table lAll were subsequently used to measure the 
relative SN 2008S magnitude for each observation. The UBVRI 
magnitudes of SN 2008S are reported with their uncertainties esti- 
mated by combining in quadrature the error of the photometric cal- 
ibration and the error in the PSF fitting in Table [5] The responses 
of the SXVF H9 camera used by A. Arbour, and the ST-9E/9XE 
camera used by M. Mobberly peak, respectively, in the V and R 
bands. Therefore, although strictly these cameras were unfiltered, 
we nevertheless list the magnitudes obtained in the V or R columns 
of TableJU 



3.2 Near Infrared data 

The near infrared (NIR) photometry was obtained with the TNG 
with NICS and with the AZT 24 telescope with SWIRCAM, in 
the JHK filters (Table [4}. The NIR images were reduced using 
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Figure 2. Finding chart for the local sequence of stars employed for the 
optical photometric calibration. The numbers are adopted from Pozzo et al. 
(2006). 



standard IRAF routines, with the jittered exposures first median- 
combined to obtain sky images in each band. Jittered images were 
then sky subtracted, registered and finally combined. The instru- 
mental magnitudes were measured on the combined images with 
SNOoPY package. Photometric calibration was carried out via rel- 
ative 2MASS photometry of the same local sequence stars used for 
optical data. The NIR magnitudes of SN 2008S are listed in Tabled 



3.3 Photometric evolution and bolometric lightcurve 

In Fig.[3](top panel) the UBVRIJHK light curves of SN 2008S 
are illustrated. These are characterised by a broad peak and a sub- 
sequent slow decline. In the R and V bands there is clearly a fast 
rise to peak shown by early observations from a range of am- 
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Table 5. Optical, Near Infrared and unfiltered photometry of SN 2008S. 



JD a 


ph° 


U 




V 


R. 


I 


J 


H 


K 


Instrument 


490.3 


4 






17 25 4- 20 


16.66zb 0.50 










MX916+ ICX424 unfiltered 


498.3 


12 






16.97 zb 0.30 












MX916 unfiltered 


501.3 


15 






16.95 zb 0.08 


16 34 _|_ o.05 










SXVF-H9+ST-9E/9XE unfiltered 


501.4 


15 








16.35 zb 0.01 










NOT 


502.2 


16 




17.92 zb 0.16 


16.87 zfc 0.12 


16.32 zb 0.08 


15.87 zb 0.07 








MSK 


503.3 


17 






16.88 zb 0.10 


16.39 zb 0.05 










SXVF-H9+ST-9E/9XE unfiltered 


505.3 


19 






16.98 zb 0.14 


16.33 zb 0.06 










SXVF-H9+ST-9E/9XE unfiltered 


506.3 


20 








16.26 zb 0.06 










SXVF-H9+ST-9E/9XE unfiltered 


507.8 


22 


17.89 zb 0.04 


I7.77 -|- o.03 


16.91 zb 0.01 


1 6 35 4 02 


15.79 zb 0.02 








TNGD 


508.3 


22 






17104014 


16.36 zb 0.06 










SXVF-H9+ST-9E/9XE unfiltered 


508.7 


23 




17.81 zb 0.30 


16.90 zb 0.50 


16 37 ± q 04 


15.86 zb 0.07 








CAHAT 


509.3 


23 






17.12 zb 0.14 


16.44 zb 0.1 1 










SXVF-H9+ST-9E/9XE unfiltered 


510.3 


24 






17.10 zb 0.14 












SXVF-H9 unfiltered 


512.3 


26 






17.08 zb 0.12 


16.42 zb 0.09 










SXVF-H9+ST-9E/9XE unfiltered 


514.3 


28 






16.99 zb 0.13 












SXVF-H9 unfiltered 


514.7 


29 




17.91 ± 0.03 


17.00 zb 0.01 


16.48 zb 0.01 


15.91 zb 0.01 








CAO 


515.3 


29 






17 06 zb 0.40 












SXVF-H9 unfiltered 


515.8 


30 


18.26 ± 0.10 


17.89 ± 0.03 


17.01 zb 0.02 


16 47 -j- o.Ol 


1591 4 01 








LT 


517.3 


32 








16.52 zb 0.08 


1 6 03 4- 1 3 








MKS 


522.7 


37 




17.98 ± 0.01 


1710 zb 0.01 


16.57 zb 0.01 


16.11 zb 0.02 








CAHAT 


523 2 


37 




18 15 ± 14 


17 08 zb 07 


16 61 ± 04 


16 20 zb 06 








MSK 


523.3 


37 






17 05 -J- 0.08 


16.54 zb 0.07 










SXVF-H9+ST-9E/9XE unfiltered 


524 7 


39 


> 18.70 


l 8 04 4 06 


17 18 4- 09 
1 / .to zn vt.vtz. 


1 67 _|_ n [-19 


16 10 4- O 01 
1 u. 1 \j zn u.u 1 








LT 


525.7 


40 






17.13 zb 0.06 












SXVF-H9 unfiltered 


526.8 


41 












15.20 zb 0.09 


15.16 zb 0.09 


14.95 zb 0.1 


TNGN 


528 7 


43 






17.22 zb 0.11 












SYVF-HQ nnfiltereH 
ja v r n; Liiuiiieicu 


530 3 


44 






17 7« _|_ n 19 












SXVF-H9 unfiltered 


530.8 


45 


18.90 zb 0.15 


18.22 zb 0.05 


17 37 _|_ o.03 


16.75 zb 0.03 


16.20 zb 0.02 








LT 


532.7 


47 




1 iq -I- n no 
10.J7 zn u.\jz 


17^8 4- n 1 


ui.oj zn u.ui 


1 £ 96 _i_ 09 








CAHAT 


533.3 


47 








[6.65 zb 0.1 1 










ST-9E/9XE unfiltered 


537.6 


52 




18.46 ±0.11 


17 46 -j- o.05 


16.91 zb 0.03 


16.45 zb 0.04 








SLV 


539.7 


54 




18.59 zb 0.01 


17.57 zb 0.02 


16.99 zb 0.02 


16.41 zb 0.01 








LT 


543.6 


57 








17.15 zb 0.05 


16.55 zb 0.05 








SLV 


544.7 


59 








17 14 _|_ 0.04 










NOT 


546.6 


61 






17.88 zb 0.14 


17.27 zb 0.10 


16.70 zb 0.07 








SLV 


546.7 


61 












15.75 ± 0.10 






AZT24 


554 2 


68 




in 9^ -)- n OS 


17 96 zb 01 


17 35 _|_ 02 


16 70 zb 01 








LT 


555 5 


69 








17 31 _|_ q 21 










SLV 


556.6 


71 






18.18 zb 0.09 


17.44 zb 0.04 


16.81 zb 0.05 








SLV 


556 6 


71 












1 6 00 +0 OQ 


1 S 69 4 O OQ 


15 39 zb 15 


AZT24 


556 7 


71 




1 q an _|_ n 00 


18 16 4- 09 


17 AT. _i_ 01 


1 6 7» -1-0 01 
lu. 1 zn u.u 1 








CAO 


557 7 


7-1 






1 8 1 ^ 4- 04 
io.u zn u.ut 


1 7 a~\ _i_ 04 
1 ; .tj zn u.u4 


1 6 77 -l-OO^ 
in. / / zn vi.kij 








LT 


570.7 


84 


20.92 zb 0.10 


20.16 zb 0.07 


18.63 zb 0.03 


17.82 zb 0.02 


17.1 1 zb 0.01 








LT 


572.7 


86 


20.93 ± 0.14 


20.17 zb 0.04 


18.78 zb 0.02 


17 07 _i_ 01 


17 i6 + 01 








LT 


573 6 


87 


70 o-i _|_ n 


70 1 8 4 OQ 


1 8 7Q 4 fll 
10./7 zn uajj 


1 7 QO -1- 01 
1 / .y\f zn u.ui 


17 14 4-0 01 
1 / . 14 zn u.u 1 


1 6 37 +0 OS 
iii.ji znu.uo 


1 c Q6 4 08 


1 S 6S 4 O OQ 


TNGD TNGN 


574.6 


88 




20.13 zb 0.07 


18.83 zb 0.03 


18.03 zb 0.02 


17.20 zb 0.02 








LT 


576.6 


90 




20.34 zb 0.14 


18.82 zb 0.04 


18.10 zb 0.03 


17.32 ± 0.01 








LT 


577.6 


91 


21.01 ± 0.30 


20.36 zb 0.01 


18.97 zb 0.03 


18.12 zb 0.02 


I7 33 _|_ q_02 








LT 


579.5 


93 






19.10 zb 0.10 


18.18 zb 0.10 


I7 33 _j_ q_07 








MSKL 


582.6 


96 






19.15 zb 0.05 












CAHAT 


585.4 


99 






19.24 zb 0.30 


18 35 4- 09 


1 7 54 4- 08 








MSKL 


590 4 


105 








10-4J zn u. 1 1 


1 7 67 -1- 1 
1 / .u / zn u. iv 








MSKL 


590.7 


105 






19.51 zb 0.01 


18.55 zb 0.01 










LT 


592.6 


107 




20 88 4 09 


19.57 zb 0.06 


18.61 zb 0.03 


17 g0 zb 0.02 








CAO 


593 6 


108 


91 68 o 90 

z 1 .uo zn u.zu 


7 1 07 4 1 3 


19 65 zb 05 


1 g 67 zb 02 


1 7 07 -1- O 09 








TNGD 


601 5 


1 1 5 








19 10 zb 2 










WOT 


602.2 


1 16 






20.06 zb 0.09 


19.15 zb 0.04 










LT 


602.5 


1 17 








19.17 zb 0.47 










MSKL 


604.5 


118 








19.28 zb 0.37 


18.10 zb 0.26 








CRM 


607.5 


122 






20.26 zb 0.26 












CRM 


607.7 


122 






20 27 4- 06 


19.26 zb 0.04 










LT 


608 1 


123 






70 7« _|_ n 1 1 


in -is. _|_ 0^ 










LT 


613.4 


127 








19.39 zb 0.19 


18.45 zb 0.15 








MSKL 


618.5 


132 






20.58 zb 0.21 












SAO 


619.5 


133 








19.45 zb 0.20 










WOT 


621.5 


135 






20.65 ± 0.27 


19 49 _|_ o.20 


18.55 zb 0.08 








CAO 


627.6 


142 






20.93 zb 0.13 


19.67 zb 0.07 










LT 


629.6 


144 












17.68 zb0.06 


16.72 zb 0.05 


16.12 zb 0.05 


TNGN 


637 4 


1 5 1 






70 Q4 -1- 14 
i^j.74 zn 14 


IQS7 + nfl(i 
1 3.0 / zn uaju 


1 8 QS 4- 0^ 
10.7J zn u.uj 








LT 


642.6 


157 






21.01 zb 0.21 


19.94 zb 0.07 


1 9 05 4- 04 

t/'UJ _L U.U4 








LT 


646.6 


159 




22.47 zb 0.22 


21.08 zb 0.14 


20.05 zb 0.09 










CAHAT 


648.4 


162 






21.18 zb 0.30 


20.02 zb 0.07 


19.07 zb 0.07 








LT 


651.5 


165 






2i 30 4j 0.30 


20.06 zb 0.05 


19 07 -J- 0.05 








LT 


658.5 


172 






21.34 zb 0.28 


20.15 =b 0.08 


19.13 =b 0.05 








LT 


659.5 


174 












17.77 =b0.06 


16.75 =b 0.08 


16.05 =b 0.08 


TNGN 


667.5 


181 








20.21 zb 0.26 


19.24 zb 0.32 








LT 


672.4 


186 






21.52 =b 0.15 


20.27 =b 0.05 


19.34 =b 0.05 








LT 


677.4 


191 






21.53 =b 0.23 


20.28 =b 0.06 


19.33 =b 0.04 








LT 


697.5 


212 












17.99 =b0.10 


16.78 =b 0.09 


16.00 =b 0.08 


TNGN 


714.4 


228 






22.09 =b 0.30 


20.75 =b 0.11 


19.65 =b 0.10 


18.14 =b0.10 


16.89 =b 0.1 


15.90 =b 0.1 


TNGN,TNGD 


782.4 


296 






22.85 =b0.30 


21.30=b 0.20 


20.30 =b 0.10 








LT 


790.4 


304 












18.85 =b 0.10 


17.20 =b0.05 


15.82 =b0.05 


TNGN 
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Table 6. Epochs, apparent and absolute magnitudes of light curve maximum 
in the BVRI bands. 



Filter 


JD a 


ph b 




M max 


B 


509 ± 2 


23 


17.83 ±0.05 


-13.76 ±0.16 


V 


505 ± 2 


19 


16.95 ±0.05 


-13.97 ±0.16 


R 


503 ± 2 


17 


16.33 ±0.05 


-14.17 ±0.16 


1 


502 ±3 


16 


15.85 ±0.05 


-14.20 ±0.16 



a JD - 2454 000.00 

6 Phase is in days after the explosion epoch JD 2 454 486 ± 4. 



Table 7. Decline rates in the BVRIJHK bands. 



Filter 


phi" 


7l(mag/100d) 


ph 2 a 


72 (mag/1 OOd) 


B 


50-100 


4.5 ±0.10 


100-160 


2.8 ±0.10 


V 


50-120 


4.0 ±0.05 


140-300 


1.3 ±0.06 


R 


60-120 


3.4 ±0.05 


140-300 


1.0 ±0.05 


1 


60-120 


2.8 ±0.05 


140-300 


0.9 ±0.06 


J 


40-120 


2.4 ±0.10 


140-310 


0.7 ±0.05 


H 


40-120 


1.5 ±0.10 


140-310 


0.4 ±0.10 


K 


40-120 


1.5 ±0.09 


140-310 


-0.2 ±0.06 



a Phase is in days after the explosion epoch JD 2 454 486 ± 4. 



ateur telescopes. The peak occurs progressively earlier from the 
blue to the red bands (Table [6j. The absolute magnitudes at max- 
imum (Table [6j, calculated adopting /j, = 28.78 ± 0.08 mag and 
correcting for Galactic (Av = 1.13 mag) and internal extinction 
(A v = 1 mag), reveal that SN 2008S is brighter by 1-1.5 mag than 
the NGC 300 OT2008-1 and M 85 OT2006-1 transients with which 
it has been compared (see Sect. [73}. 

The BVRI light curves show very similar temporal evolu- 
tions with three phases characterized by a different decline rate: a 
broad peak (about two weeks), a phase of steeper decline (71 in 
Table |7J starting about 60 days after explosion, more pronounced 
at shorter wavelengths, and a flattening (72) after 140 days, more 
evident at longer wavelengths. Due to the faintness of SN 2008S 
we do not have U and B band observations at the late phases to 
obtain an accurate estimate of the decline rate. In the NIR bands 
the decline rates are very slow until 120 days and show further a 
flattening (or a slight increase in K band) after this epoch. 

The time evolution of the B - V, V - R, V - I, J - H and 
J — K colours of SN 2008S is shown in Fig. [3] All optical colours 
become progressively redder until about 200 days after the explo- 
sion. Subsequently the colours do not show significant evolution. 
This trend is most evident for the V — I colour. The NIR colours 
show a different evolution after about 200 days after the explosion: 
J — H is slightly increasing while J — K is steepening. 

A "bolometric" light curve (Fig. [3] bottom panel) was 
obtained by first converting UBVRIJHK magnitudes into 
monochromatic fluxes per unit wavelength, then correcting these 
fluxes for the adopted extinction (Ay = 2. 13 mag) according to 
the extinction law from lCardelli et all dl989h . and finally integrat- 
ing the resulting SED over wavelength, assuming zero flux at the 
integration limits (the blue end of the U band and the red edge of 
the K band). We estimated the flux only for the phases in which V 
band observations were available. The photometric data in the other 
bands were estimated at these phases by interpolating magnitudes 
in subsequent nights. 

During the period 140-290 days, the UBVRIJHK "bolo- 



metric" light curve tail shows a decay rate of 0.88 ±0.05 mag/lOOd , 
very similar to that of 56 Co (1.023 mag/lOOd from lHuo et alll987l) . 
while the UBVRI "bolometric" lightcurve is slightly steeper 
(1.3 ± 0.05 mag/1 OOd). Assuming that radioactive material was 
powering the late time photometric evolution of SN 2008S, we tried 
to estimate the 56 Ni mass synthesized by SN 2008S by compar- 
ing its U BVRIJHK quasi-bolometric light curve with that of SN 
1987A from 140 to 300 days after explosion and assuming similar 
7— ray deposition fraction: 

M( 56 Ni)o8S = M( 56 Ni) 87A x £251 (1) 

where M( 56 Ni) 87A is the mass of 5B Ni produced by SN 1987A, 
Loss is the luminosity of SN 2008S and Lg7A is the luminosity of 
SN 1987A (at a similar epoch) also obtained from UBVRIJHK 
data. We adopt M( 56 Ni) 87 A = 0.0 73 ± 0.012 Mp whi ch is the 
weighted mean of v alues given by lArnett & FvJ j 1989b and by 
iBouchet et all dl99lh . For SN 2008S we obtain a M( 56 Ni) of 
0.0014±0.0003 Mq where the error includes both the uncertainties 
in the assumed distance of SN 2008S and in the 56 Ni mass of SN 
1987A. If we considered only optical data for both SN 2008S and 
SN 1987A we would obtain a M( 56 Ni) of 0.0011± 0002 M w . We 
estimated the 56 Ni mass also using the method of lHamuvl (2003) 
assuming that all the 7— rays resulting from the decay of 56 Co into 
56 Fe are fully thermalized: 

M( 56 Ni)o8S = 7.866X 10' 44 Lexp[ {t - * o)/(1 + z) - Tm ]M (2) 

TCo 

where to is the explosion epoch, t Ni =8.76 days the life time of 
56 Ni and r Co = 111.26 days is the life time of 56 Co. Using this 
method we estimated M( 56 Ni) for each point of the radioactive tail 
and the average of these estimates gives an ejected 56 Ni mass of 
~ 0.0012 Mq. We remark that the estimated value of 56 Ni mass 
has to be considered as an upper limit if the ejecta-CSM interaction 
or an IR echo contributes to the tail of the bolometric light curve. 

3.4 Data at other wavelengths 

SN 2008S was serendipitously observed with Spitzer on 2008 Feb 
6.8 UT during scheduled obser vations of the nearb y SN 2002hh 
showing a strong MIR emission ( Wesson et alj2008l) . Images were 
obtained with IRAC (3.6 fim, 4.5 fim, 5.8 ^im and 8.0 /im) on 
day 17.3 and MIPS (24 /im) on day 18.0, uniquely early epochs 
for coverage of this wavelength region. The Swift satellite also 
observed SN 2008S in the optical, UV and X-ray bands on Feb. 
4.8, 6.0, and 10.5 UT and detected it only in the optical im- 
ages dSmithetal]|2008h . SN 2008S was not detected on Febru- 
ary 10.62 UT at radio frequ encies with the Very Large Array by 
Chand ra & Soderberd (2008). The contribution to the overall en- 
ergy budget of SN 2008S at about 20 days after the explosion seems 
to be negligible at UV wavelengths but it is considerable at MIR 
wavelengths. In Sect. 14. 11 we shall demonstrate that the large MIR 
luminosity must have been caused by an IR echo from pre-existing 
circumstellar dust. 

4 SPECTRAL ENERGY DISTRIBUTION EVOLUTION 
4.1 Up to day 120 

The SED evolution of SN 2008S is shown in Fig. 0] During the 
first ~ 120 days, the optical-NIR fluxes can be well reproduced by 
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Figure 3. Top panel: UBVRIJHK light curves of SN2008S. Middle panel: evolution of the B - V, V - R, V - 1, J - H, J - K colours. The magnitudes 
have been corrected for the adopted extinction according to the Cardelli extinction law. Bottom panel: the U BVRIJHK and UBVRI quasi-bolometric light 
curves. The 56 Co decay slope and the slope of UBVRIJHK and UBVRI curves are shown. Phase is in days after M 8$Q>&&?t$s8P$§ 
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Table 8. Parameters for a single blackbody fit (until 108 days after explosion) and for a two-component (hot and warm blackbodies) fit to the UBVRIJHK 
fluxes of SN 2008S. 



pn 


T, 


K hot 


T , 


T 

1 warm 


D 

R warm 


T 

^warm 


T f\ 

L warm' L'tot 




\Z 
IV 


i ril4 f.™ 

iu cm 


iu eig s 


tv 


1 nl4 r-m 

iu cm 


i n40 arn c — 1 

iu erg s 


per cent 


17 


8380 ± 150 


1.93 ± 0.08 


13 ±0.2 










22 


8400 ± 120 


1.90 ±0.08 


12 ±0.1 










39 


7600 ± 100 


1.89 ±0.07 


8.4 ± 0.1 










71 


6300 ± 90 


1.77 ±0.08 


3.5 ±0.4 










108 


5220 ± 70 


1.50 ±0.07 


1.2 ±0.1 










122 


5000 ± 500 


1.30 ±0.20 


0.7 ±0.3 


1400 ± 500 


5.0 ±4.0 


0.08 ± 0.1 


10 


142 


4920 ± 160 


1.09 ± 0.08 


0.5 ±0.1 


1430 ± 200 


5.3 ±0.5 


0.11 ± 0.1 


25 


172 


4670 ± 180 


0.99 ± 0.10 


0.3 ±0.08 


1480 ± 150 


6.3 ±1.0 


0.14 ±0.09 


32 


228 


4430 ±170 


0.85 ± 0.08 


0.2 ±0.05 


1413 ± 100 


8.4 ±2.0 


0.20 ±0.1 


50 


296 


4380 ± 160 


0.64 ±0.06 


0.1 ±0.03 


1200 ± 60 


13 ±0.2 


0.28 ± 0.08 


74 



1 Phase is in days after the explosion epoch JD 2 454 486 ± 4. 
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Figure 4. Temporal evolution of the observed SED of SN 2008S. Lines 
show the blackbody fits to the SED. Observed fluxes are corrected for 
the adopted extinction (Ay = 2. 13 mag) according to the extinction law 
from Cardelli et al (1989). Phase is in days after the explosion epoch (JD 
2 454 486). 



a single hot blackbody. The blackbody temperatures, radii and lu- 
minosities for some epochs are shown in Table [8] and plotted in 
Fig.H] 

The temperature fell monotonically from ~ 8300 K to 
~ 5000 K during this time. The radius declined slowly from about 
1.9 x 10 14 cm to 1.3 x 10 14 cm, suggesting that the blackbody 
surface was defined by the photosphere through which the ejecta 
flowed. 

Remarkably, a strong MIR excess was observed at ~ 18 days 
after explosion by IWesson et al. | ilOOA The SED obtained with 
the Spitzer data was fitted by a 500 K blackbody modified by a dust 
emissivity that is inversely proportional to the wavelength, giving a 
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Figure 5. Temporal evolution of the parameters of the hot (left) and warm 
(right) blackbodies fitting the SED of SN 2008S. Phase is in days after the 
explosion epoch (JD 2 454 486). 



luminosity of 21 x 10 Lq. We have re-measured the Spitzer im- 
ages, using the post-basic calibrated data (PBCD) products pro- 
vided by the Spitzer pipeline. The 8 tim and 24 fim images showed 
strong irregular background at or near the location of the SN. We 
therefore subtracted serendipitously obtained pre-explosion "tem- 
plates" from the post-explosion images before proceeding with our 
flux measurements. The IRAC templates were taken from Spitzer 
programme 30292 (RI. Meikle) and the MIPS template from pro- 
gramme 0230 (P.I. Barlow). The image matching and subtraction 
was performed as implemented in the ISIS v2.2 image subtrac- 
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Figure 6. Fit of two blackbodies to the SN 2008S SED from optical to 
MIR at 17.3 days after explosion. The NIR data are obtained extrapolating 
data from phase 41. Observed fluxes are corrected for the adopted extinc- 
tion (Ay = 2.13 mag) according to the extinction law from Cardelli et al 
(1989). 



tion package dAlardll2000l), and modifi ed in a manner analogous 
to that described in iMeikle et alJd2006l) . Aperture photometry was 
performed on the background- subt racted IRAC and MIPS images 
using the Starlink package GAIA torapeill200(]|) . A circular aper- 
ture of radius 5" was used for the photometry. This was chosen as 
a compromise between minimising the effects of the residual irreg- 
ular background emission at the SN location and minimising the 
the size of aperture correction needed in the final flux determina- 
tion. The aperture radius corresponds to a distance of ~ 140 pc at 
SN 2008S. Residual background in the template-subtracted IRAC 
and MIPS images was measured and subtracted by using a clipped 
mean sky estimator, and a concentric sky annulus having inner and 
outer radii of 1.5 and 2 times the aperture radius, respectively. Aper- 
ture corrections were derived from the IRAC and MIPS point re- 
sponse function frames available from the Spitzer Science Cen- 
ter, and ranged from xl.04 at 3.6 /im to x2.12 at 24pm. For 
IRAC, the aperture was centred by centroiding on the SN image. 
For MIPS, the aperture was centred visually on the SN, checking 
against the WCS coordinates. We found fluxes of 1.60 ± 0.02 mjy 
at 3.6 ^m, 1.97 ± 0.03 mjy at 4.5 ^m, 3.06 ± 0.07 mjy at 5.8 
and 4.15 ± 0.05 mjy at 8 ^m. The errors are statistical only and 
do not include calibration uncertainties which may amount to an 
additional ~ ±5%. The image subtraction at 24 /im left a substan- 
tial residual. We therefore also carried out aperture photometry on 
the unsubtracted 24 image. Consistent results were obtained for 
the subtraction and non-subtraction procedures. The mean flux de- 
rived from the two methods is 0.7 ± 0.1 mjy Our IRAC and MIP S 
fluxes are consistent with those reported by [Wesson et alj d2008l) . 
However, their errors are at least several times larger. 

To investigate the nature of the MIR source, we carried out a 
simultaneous fit of two blackbodies to the SED at 17.3 days, and 
this is shown in Fig. [6] It can be seen that the blackbodies provide a 
plausible representation of the data. The hot blackbody has a tem- 
perature of 8076 ± 150 K and a radius of (2.1 ±0.1) x 10 14 cm, 
corresponding to an expansion velocity of 1430 kms -1 . These es- 
timates are consistent with emission from the hot photosphere. The 
warm blackbody has a temperature of 585 ± 5 K. This temperature 
together with the reasonable match to the blackbody leads us to pro- 
pose that the MIR emission must be due to warm dust. The radius 
of the warm blackbody is (9.9 ±0.4) x 10 15 cm. For the SN to pro- 



duce a surface of this radius in just 17.3 days would take a velocity 
of 66, 200±2700 kms" 1 . Such an enormous velocity immediately 
rules out newly-formed ejecta dust as the source of the MIR emis- 
sion. The high velocity also rules out collision of the fastest moving 
ejecta with a dusty CSM. The only viable alternative appears to be 
an IR echo from CSM dust. To explore this possibility we have 
matched an IR echo model to th e MIR SED. Details of the model 
are given in Mei kle et al This assumes a spherically sym- 

metric cloud of carbon grains centred on the SN, with a concentric 
dust-free cavity at the centre. For simplicity, a single grain radius, 
a, is adopted. For ease of computation, we assumed that the grain 
material was amorphous carbon where, for wavelengths longer than 
2™, the grain absorptivity/e missivity can be well-app roximated as 
being proportional to A -1 15 dRouleau & Martini 1 99 lh . For shorter 
wavelengths, an absorptivity/e missivity of unity was use d. The ma- 
terial density is 1.85 gem -3 dRouleau & Martin 1 199 lh . Free pa- 
rameters are the grain size, grain number density, the CSM radial 
density law, the CSM extent and the size of the concentric dust- 
free cavity. The outer limit of the CSM was set at 10 times that 
of the cavity radius, although this parameter is not critical. The in- 
put luminosity is a parametrized description of the UBVRIJHK 
"bolometric" light curve (Fig. |3). Given the apparently low tem- 
perature of the dust, it is reasonable to assume that the IR echo 
made only a small contribution to the NIR flux implying a negligi- 
ble overestimate of the bolometric light curve. 

It was found that a range of parameters were able plausibly 
to reproduce the MIR flux. Matches were found for grain radii of 
0.001-0.5 fim. The corresponding cavity range was 1 x 10 17 cm 
to 3 x 10 cm and the dust mass range was 10" 2 - 10" 3 M Q . 
Such masses of dust are entirely plausible within the CSM of a 
red (super)giant. For a typical wind velocity of 20 kms _1 and a 
dust-gas mass ratio of 0.01, the CSM dust masses we require in 
the IR echo model would be produced by a mass loss rate of a few 
times 10 Mq yr _1 . All the progenitor circumstellar dust within 
a particular radius would have been evaporated in the SN explo- 
sion. Given the peak luminosity of the SN, the evaporation radius 
for carbon grains is 3.5 x 10 15 cm. Silicate grains would allow 
an evaporation cavity several times larger. However, there is lit- 
tle sign of an excess in the 8 /im band which we would expect 
from optically thin silicate grains. Nevertheless, even to account 
for a cavity radius of 3 x 10 16 cm we must invoke some sort of 
episodic mass-loss history. The echo model spectrum is shown in 
Fig. [7] left-hand panel. Also shown are the hot blackbody spectrum 
obtained in the two-blackbody fits described above, and the com- 
bined hot blackbody and IR echo spectrum. This is compared with 
the optical-NIR-MIR observations at this epoch. The earliest NIR 
data were obtained at 40 days and so the 17.3 days fluxes were 
estimated by extrapolation of the NIR light curves. The extrapola- 
tion was performed by assuming constant V — J, V — H, V — K 
colours during the first month. (If, instead, we assume a linear evo- 
lution of these colours, then the fluxes in the NIR bands would be 
smaller by 6%, 22% and 27% in J, H, K respecti vely.) The hot 
blackb ody was reddened using the extinction law of ICardelli et al.1 
( 1989) with Ay = 2.13 mag. The total optical depth in the opti- 
cal region through the CSM in the IR echo model is about 0.20, 
or Ay — 0.22 mag. The reddening effects of the CSM dust are 
included in the model. Consequently the final IR echo model was 
reddened with Ay — 2.13 — 0.22 = 1.91 mag. Inspection of Fig. [7j 
shows that a satisfactory match to the data is obtained. 

While the echo model easily reproduced the MIR flux, it also 
produced a certain level of NIR emission, with K being the bright- 
est of the observed NIR bands. This, therefore, provides an addi- 
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tional constraint on the echo model i.e. the model has to be consis- 
tent with the NIR flux and its evolution. As explained above, satis- 
factory fits to the early-time optical-NIR SEDs were obtained using 
a single blackbody. There was little evidence of a NIR excess up to 
around 120 days. However, the blackbody fits are subject to er- 
ror. We therefore investigated whether or not the predicted K band 
fluxes of the echo model could be consistent with the uncertainties 
in the observed values and in the hot blackbody fit. For a number 
of epochs (17.3, 38.7,71.7, 107.6, 116.2, and 151.4 days) we com- 
pared the echo model with the residual K band fluxes derived by 
subtracting the hot blackbody flux to the observed one. The first five 
epochs were all consistent with a zero residual, but with uncertain- 
ties of ~ ±0.1 mjy. We then selected the echo model which yielded 
the smallest K band fluxes and still produced a satisfactory match 
to the MIR fluxes. This was achieved with a 0.5 /im grain radius and 
a dust-free cavity of 3 x 10 cm. The dust mass is 1.4 x 10" 3 M o . 
The right hand panel of Fig. [7] shows the residual fluxes at 2.2 /im 
together with the 2.2 /im light curve from the echo model (coloured 
blue). The error bars on the residual points were derived as a com- 
bination of the uncertainty in the observed/extrapolated points and 
the uncertainties in the hot blackbody fits. It can be seen that the 
model flux never exceeds the residual by much more than 1 a. The 
marginal negative shift of the first three residual points may be due 
to slightly imperfect blackbody fits such as might be caused by the 
actual extinction law being different from the one adopted. We de- 
duce that it is quite plausible that the K band flux from an IR echo 
would have been undetected, hence accounting for the near-zero 
residuals. The later three residual points show the gradual emer- 
gence of the NIR excess, discussed below. We conclude that the IR 
echo model accounts for the MIR flux as well as being consistent 
with all the shorter wavelength fluxes. 

Also shown in Fig.|7](right hand panel) is the 8 fim light curve 
predicted by the IR echo model. We note that neither the 2.2 /im 
nor 8 /im light curves exhibit the flat top or plateau characteristic 
of the IR echo phenomenon. This is because, owing to the excep- 
tionally low bolometric luminosity of SN 2008S, the radius of the 
dust free cavity had to be small in order to attain a sufficiently high 
dust temperature. Indeed, the radius of the cavity in the model il- 
lustrated here is only 12 light days (3 x 10 16 cm). Consequently the 
decline of the light curves are actually dominated by the character- 
istic decline timescale of the input bolometric light curve, which is 
about 70 days per e-fold, together with the r ~ 2 dust density decline. 

The final MIR observation of SN 2008S was at about 
180 days, but at the time of writing, the data are not yet out of 
embargo (due 2009 August). Nevertheless, we can use the IR echo 
model to predict the fluxes at this epoch. We find that the 8 /im 
flux at 180 days is 1.3 mJy for the model discussed above, rising to 
around 5 mJy for models having a grain radius ~ 0.01 /im. Thus, 
the observed 1 80 day MIR fluxes should provide constraints for the 
grain size. 



4.2 After day 120: dust formation? 

To fit the SED at later epochs we need the sum of a hot and a 
warm blackbody. We find statistical evidence for the emergence of 
a warm component as early as 120 days, viz. ABIC = BIChot — 
BIChot+warm — 1.5. This becomes stronger with time, reaching 
ABIC ~ 5.7 by 135 days. The hot component was still dominant 
at 172 days but by 296 days it contributed only 26% of the total lu- 
minosity. The temperature of the warm component stayed roughly 
constant at 1400-1500 K until at least 228 days, but by 296 days 
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Wavelength (/im) Epoch (days) 

Figure 7. The left hand panel shows the day 17.3 hot blackbody spec- 
tram (green), the IR echo spectrum (red) and the combined spectrum (blue). 
These are compared with the day 17.3 optical-NIR points and the day 17.3 
MIR points. The hot blackbody has a temperature of 8076 ± 150 K and 
a radius of (2.1 ± 0.1) X 10 14 cm, reddened using the extinction law 
of Cardelli et al. (1989) with Ay = 2.13 mag. The total optical depth 
through the CSM in the optical region is about 0.20, or Ay = 0.22 mag. 
The reddening effects of the CSM dust are included in the model. Conse- 
quently the final IR echo model was reddened with Ay = 2.13 — 0.22 = 
1.91 mag. For the IR echo model, a single grain size was used with radius 
a = 0.5 /tm. The grain material is amorphous carbon where, for wave- 
lengths lon ger than 2ir/a, the grain absorptivity/emissivity is proportional 
to A -115 iRouleau & Martin 1991). For shorter wavelengths, an absorp- 
tivity/emissivity of unity was used. The material density is 1.85 gem - 3 
iRouleau & Martini 1 199 lb . The CSM extended from 3 x 10 16 cm to 3 X 
10 17 cm, with the dust density declining as r~ 2 (steady wind). The dust 
mass is 1.4 X 10 -3 Mq. The right hand panel shows the residual fluxes at 
2.2 fim (blue points) at a number of epochs together with the 2.2 /tm light 
curve from the IR echo model (blue). For the first five points it can be seen 
that the model flux never exceeds the residual by much more than lcr. The 
latest three residual points show the gradual emergence of the NIR excess 
(see text). Also shown is the IR echo light curve at 8 /tm (red). 

had cooled to 1200 K. The luminosity and the radius monotonically 
increased during the 122-296 days period (Fig. [5}. 

A clear NIR excess was also observed at late phases in several 
other SNe(1997ab, 1979C, 1980K, 1982E, 1982L, 1 982R, 1985L 
1993 J , 1994Y, 1995N, 1998S, 2005ip, 2006jc) (e.g. lGerardv eTail 
2002). The simplest explanation for the IR excess in both SN 2008S 
and these other SNe is thermal emission from warm dust heated by 
the SN. However, different locations, origins and heating mech- 
anisms are possible. Specifically, IR excesses may be due to (a) 
newly-formed dust in the ejecta, heated by radioactive decay or a 
reverse shock, (b) newly-formed dust in a cool dense shell (CDS) 
formed by the SN shock/CSM collision and heated directly by the 
shock, by absorption of X-rays from the interaction region or by the 
bolometric light curve of the SN (a type of IR echo), (c) ejecta col- 
lision with pre-existing circumstellar dust or (d) an IR echo from 
pre-existing circumstellar dust heated by the SN bolometric light 
curve. These scenarios do not necssarily exclude one another. For 
example, in the case of SN 2006jc the NIR excess was explained 
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bv lMattilaetaD ( l2008h as a combination of IR echoes from newly- 
formed CDS dust and from pre-existing dust in the CSM. 

In the case of SN 2008S, we can immediately rule out an 
IR echo, whether from CDS or CSM dust, as being the cause of 
the late-time NIR excess. The bolometric light curve is too weak 
by a factor of ~ 100 to account for the high dust temperatures and 
luminosities seen during the 200-300 days period. Moreover, for 
an IR echo scenario we would expect the phase and rise-time of the 
echo flux to be determined by the temporal behaviour of the bolo- 
metric light curve, which is about 35 days per e-fold beginning at 
days (cf. Fig. [7] right-hand panel). Yet the NIR excess does not 
appear until at least 120 days and has a characteristic risetime of 
about 100 days. 

This leaves dust heated by radioactivity or by energy from 
shock interaction with CSM. We first estimated the mass of such 
dust required to reproduce the NIR flux. We employed the es- 
cape probability model described in Mei kle et al.l d2007h . Setting 
the expansion velocity of the dust cloud at 1000 ± 50 tons" , 
we found that for silicate dust grains in the ejecta the NIR ex- 
cess was reproduced with a dust mass rising monotonically from 
about 0.2 x 10~ 5 M Q on day 122 to 1.2 x 10" 5 M Q on day 296. 
During this time, the dust temperature fell from about 1300 K to 
1 100 K and the optical depth to the centre at 2.2 /im rose from 0.2 
to 0.4. Thus, optically thin dust lying well within the velocity limit 
of refractory elements (say up to ~ 2000 kms -1 ) was able to re- 
produce the NIR excess. It follows that the derived masses indicate 
the minimum dust mass required. Reducing the size of the dust 
cloud caused it to quickly become optically thick, with the black- 
body limit being approached at a radius corresponding to about 
500 kms -1 . Clearly, in this optically thick case the NIR obser- 
vations do not restrict the mass of dust that could be concealed 
in the ejecta. If instead of ejecta dust, we assume the dust lay in 
the cloud which gave rise to the early time MIR echo, and was 
heated by ejecta collision, then the mass of dust required to repro- 
duce the NIR flux remains much the same, but the optical depth 
falls to around 0.001. Finally, for optically thin dust, replacing sili- 
cate grains with amorphous carbon grains reduces the derived dust 
mass by a factor of about 10, due to the higher emissivity of carbon 
in the NIR region. 

We conclude that a very modest mass of optically thin dust can 
account for the NIR excess flux. Moreover, it appears that the dust 
mass involved grew quite rapidly between 122 and 296 days. This 
dust may have condensed in the ejecta, heated by radioactive decay 
or a reverse shock. Alternatively, its location may have been much 
farther out in the CSM. To account for the early-time MIR flux we 
had to invoke a "dust-free cavity" of radius at least 3 x 10 16 cm, 
with a dust mass beyond this limit of at least 10 -3 Mq. However, 
to account for the later NIR excess flux, we have shown that the 
required dust mass would be barely 1% of this. Thus, it is quite 
feasible that sufficient dust to account for the NIR flux lay within 
the MIR echo cavity, with an inner limit fixed only by a dust evapo- 
ration radius which could be as small as ~ 3 x 10 cm (see above). 
For the ejecta to reach this distance in, say, 150 days, would take 
a velocity of just 2300 kms -1 . The apparently increasing mass of 
dust might have been due to dust newly forming within a CDS. Al- 
ternatively, it simply may be that the shock swept up an increasing 
mass of pre-existing CSM dust. 

Is it plausible that the NIR excess was due to new, 
radioactively-heated ejecta dust? We have already shown that this 
can be achieved with a small mass of optically-thin dust lying com- 
fortably within the limits of refractory elements. However, exami- 
nation of the energy budget reveals problems. In subsection |7.4| we 



estimate that SN 2008S produced only about x0.02 of the mass of 
56 Ni that was released by SN 1987 A. If we scale t he radioactive en- 
ergy deposition law for SN 1987 A dli et aljl993h to SN 2008S we 
find that the deposited radioactive energy exceeded that of the dust 
luminosity by x8 on day 122 falling to x 1.25 by day 228 i.e. there 
was sufficient radioactive energy to power the dust luminosity up 
to about day 228. However, by day 296, even if we assume the en- 
tire radioactive luminosity was deposited in the ejecta, it could only 
supply half of the dust luminosity i.e. for radioactivity to account 
for the warm component flux at 296 days we would have to double 
the derived mass of 56 Ni. Another possible difficulty is that, even if 
we focus on just the period 122-228 days, the fraction of radioac- 
tive energy deposited in the dust increases from 12 ~ % to 80 ~ % 
and yet the average number density of ejecta dust stayed approxi- 
mately constant. Why then did the dust-deposition fraction show 
such a large increase? It may be related to the detailed distribution 
of the grain growth and radioactive material. But, given the day 296 
energy deficit, perhaps a more plausible explanation is that at least 
part of the NIR excess had another cause. This might be reverse- 
shock heating of new ejecta dust. Alternatively it might be IR ra- 
diation from shock-heated newly-formed CDS dust or pre-existing 
circumstellar grains. We conclude that, while the earlier NIR ex- 
cess could have been entirely due to new, radioactively-heated dust 
in the ejecta, as time went by, an increasing proportion of the flux 
must have been powered by other mechanisms. 

NIR/MIR-based evidence for newly-formed dust in the 
ejecta or CDS of SNe has been reported in the cases 
of SN 1987A dMoselev et al.l 1 19891; Isuntzeff & Bouchet fl990l 



Meikle et alj Il993l : iRoche et al.l Il993h), SN 1998S jPozzo et al 
2004b . SN 2004gd dSugerman et al.l 2006; M eikle et alj 1200' 



SN 2004et footak et .3120091) and SN 2006ic JSmith i et al.ll2008l : 
|Pi Carlo et al. 2008; Mattil a et alj|2008l) . Dust masses of typically 
10 -4 -10 -3 M0 are directly observed, although larger masses may 
exist in optically thick clumps. Dust condensation in the ejecta or 
in a CDS can be also be demonstrated via the effects of the dust 
on optical radiation. In particular, it can attenuate the red wings 
of sp ectral lines, causing an apparent blue shift of the line pro- 
files fcucv et al.ll99ll;|Panziger et al.ll99ll ; lElmhamdi et al.l2003l ; 
IPozzo et alj|2004l) . However, the spectra of SN 2008S are of insuf- 
ficient resolution or wavelength precision to enable us to make this 
test. 



5 SPECTROSCOPIC DATA REDUCTION AND 
ANALYSIS 

We spectroscopically monitored SN 2008S for eight months with 
several telescopes and details of the spectroscopic observations are 
reported in Table [9] Spectra were reduced (trimmed, overscan and 
bias corrected, flat-fielded) using standard routines within IRAF. 
An optimal variance weighted extraction of the spectra was carried 
out using the IRAF routine APALL. Wavelength calibration was 
performed using the spectra of comparison lamps acquired with 
the same instrumental configuration as the SN observation. Flux 
calibration was done using spectrophotometric standard stars ob- 
served with the same instrumental set-up. Approximate spectral 
resolutions were estimated from the full-width-at-half-maximun 
(FWHM) of the night sky lines. The wavelength calibration is ac- 
curate to ±1 A for ISIS spectra, ±2 A for DOLORES and DBSP 
spectra and ±3 A for CAFOS, ALFOSC and AFOSC spectra. At- 
mospheric extinction corrections were applied using tabulated ex- 
tinction coefficients for each telescope site. The spectra of standard 
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Table 9. Journal of spectroscopic observations of SN 2008S. 



JD a ph 6 Telescope Grism or Grating Range Resolution exp. time 

A (FWHM) A s 



500.8 


15 


WHT 


R1200R 


6200-6880 


0.7 


1800 


501.3 


16 


NOT 


Gr4 


3600-8700 


19 


800 


504.3 


19 


WHT 


R158R 


5500-10000 


10 


800 


507.7 


22 


TNGD 


LR-R 


5150-10225 


11 


1800 


513.7 


28 


CAO 


Gr4 


3900-7800 


24 


3600 


514.8 


29 


WHT 


R1200R 


7025-7560 


0.8 


900 


514.8 


29 


WHT 


R300B 


3170-5350 


3.5 


900 


523.7 


38 


CAHAT 


G-200 


4950-10260 


13 


2700 


529.8 


44 


INT 


R300V 


3750-9340 


5 


1000 


532.7 


47 


CAHAT 


G-200 


4950-9750 


14 


3000 


539.8 


54 


TNGD 


VHR-V 


4650-6600 


5 


2400 


552.8 


67 


CAO 


Gr4 


3790-7790 


24 


3600 


553.7 


68 


CAHAT 


r-100 


5800-9600 


9 


3600 


JJ /.o 


1 1 


LAU 


KjX 4 


i ^ nn q a ^ n 
JjUU— o4jU 


Z4 


jOUU 


561.7 


76 


TNGD 


VHR-R 


7110-7560 


4.1 


2700 


582.6 


97 


CAHAT 


r-200 


6178-10600 


11 


3600 


588.7 


103 


CAHAT 


b-200 


3900-8850 


11 


1700 


606.7 


121 


WHT 


R316 


5520-8260 


6 


2700 


609.7 


124 


WHT 


R316 


7190-10400 


4.7 


1800 


667.7 


182 


WHT 


R158 


5400-9900 


10 


2700 


668.7 


183 


WHT 


R158 


5400-10200 


6 


2700 


741.7 


256 


HP200 


300 


3300-5600 


8 


1800 


741.7 


256 


HP200 


158 


5800-10300 


11 


1800 



a JD - 2454000.00 
b Phase in days with respect to the explosion date JD 2 454 486 ± 4. 
WHT = the 4.2 m William Herschel Telescope (WHT) with the Intermediate dispersion Spectrograph and Imaging System (ISIS); NOT= the NOT with 
ALFOSC; TNGD = the TNG with DOLORES; CAO = the Copernico telescope with AFOSC; CAHAT= the 2.2 m telescope at CAHA with CAFOS; INT= 
the 2.5 m Isaac Newton Telescope (INT) with the Intermediate Dispersion Spectrograph (IDS); HP200= the 5.08 m Hale Telescope with the Double 

Spectrograph Specs (DBSP) at Palomar Observatory. 

stars have also been used to identify telluric features and to remove 
these from the SN spectra. Spectra of similar quality obtained dur- 
ing the same night were combined to increase the signal to noise 
ratio (S/N) ratio. To check the flux calibration, BVRI magnitudes 
were estimated by integrating the spectra convolved with standard 
filter functions using the task CALCPHOT within the IRAF pack- 
age STSDAS. The spectro-photometric magnitudes were compared 
to the photometric observations and, whenever necessary, a scale 
factor was applied to match the photometric observations. The flux 
calibration is accurate to within approximately 10%. 

5.1 Line identification 

We identified Ha, H/3, H5, [Ca II] doublet (AA7292,7324) and Ca 
II triplet (AAA8498,8542,8662), Fe II (multiplets 27, 28, 37, 38, 40, 
42, 46, 48, 49, 72, 73, 74, 92, 186, 199), [O I] and Mn II (multi- 
plet 4) in emission and Na I doublet (AA5890,5896) in absorption 
in early spectra (Fig. [8] and Table [A2t . In the more recent spectra 
O I (multiplet 4, A8446) and Na I D are visible in emission, while 
the O I triplet in absorption seems to disappear. These emission 
lines are likely produced in different regions with different ioniza- 
tion states and velocities. The low ionization is likely due to high 
column density material with high optical depth. 

5.2 Spectroscopic evolution 

The spectra of SN 2008S do not show significant evolution during 
the temporal interval from 15 days to eight months after the explo- 
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sion (Fig. 1 10b, the continuum becoming progressively redder and 
fainter with time. All spectra consist of a nearly featureless contin- 
uum with superimposed strong Balmer emission lines, [Ca II] dou- 
blet and Ca II near-infrared triplet. The most remarkable change 
in the latest spectra is the presence of Na I D and O I (A8446) in 
emission. We do not detect high velocity absorption lines at any 
phase, suggesting that the ejecta are not directly visible. The lines 
do not show a P-Cygni profile typical of a SN explosion at very 
early phase but exhibit two different kinds of profiles. [Ca II] and 
Fe II lines show only a narrow asymmetric component with a red 
wing (Fig.[9]and Fig. [14). The Ha, H/3 and Ca II NIR triplet lines 
show evidence for a multicomponent profile (Fig.(9] top and middle 
panels) We focused our analysis on Ha, the [Ca II] doublet and the 
Ca II triplet since they are the most prominent lines and visible in 
all spectra estimating line parameters with the IRAF task SPLOT. 
The velocities are FWHM and those of Ha and Ca II (A8662) have 
been measured by deblending the multiple-component profile and 
assuming a Gaussian profile for each component. To measure the 
peak position at different phases we selected the spectra acquired 
with a slit width smaller than the seeing, hence considering only 
spectra with accurate wavelength calibrations. 

Hydrogen 

The intensity of Ha decreases by a factor of approximately ten 
from the earliest spectrum to the latest (Table [Tol. It follows the 
flux in R band (see Fig. lilt during the first 60 days, then shows a 
flattening (the R band light curve changes the decline slope shortly 
after). The Ha profile in the first high resolution ISIS spectrum ex- 
hibits three different kinetic components as shown in Fig. [9] The 
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Figure 8. Spectral lines observed in SN 2008S spectra. From the top to the bottom panel: spectra obtained at WHT (phase 29), at TNG (phase 22), at WHT 
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Figure 10. Time sequence of SN 2008S spectra. The spectra are corrected for Galactic extinction. The phase reported to the right of each spectrum is relative 
to the explosion epoch (JD 2 454 486) and wavelength is in the observer frame. The © symbols mark the positions of the most important telluric absorptions. 



narrow, intermediate and broad components correspond to veloc- 
ity widths (FWHM) of v n ~ 250 kms" 1 , vi ~ 1000 kms _1 and 
Vb ~ 3000 kms -1 , respectively, in this early spectrum. The Ha 
intensity and the broad component velocity width, Vb, show a simi- 
lar temporal behaviour, decreasing quickly until about 60 days after 
explosion and more slowly after this epoch (Fig.QT]and Tablelllt. 
The intermediate component velocity width, vj, shows a less abrupt 
decline. The narrow component, v n , is not resolved in several spec- 
tra at intermediate phases but appears to stay constant until 60 days 
after the explosion and is not visible after this epoch (Fig. II Hand 
Tablell It. This component of Ha seems also to be asymmetric, but 



the resolution of our spectra is not adequate to analyse its profile in 
detail. 

The Ha profile could be interpreted as broad underlying emis- 
sion due to the ejecta, an intermediate component resulting by 
shocked material behind the interaction front between ejecta and 
CSM and a narrow component from unshocked CSM. However, 
the lack of any high velocity absorption component, or P-Cygni 
profile, is a fairly strong argument that the broad emission is not 
a direct measurement of the ejecta velocities. This profile may be 
also interpreted as a narrow core with broad wings resulting from 
multiple scattering events with thermal electrons, an escape mech- 



© 0000 RAS, MNRAS 000. 000-000 



16 M.T. Botticella et al. 



Table 10. Evolution of the intensity of Ha, Ca II (A8662) and [Ca II] doublet. Intensities are in units of 10 14 erg cm 2 s 1 . 



JD a 


ph 6 


I (Ha) 


I (Ca II A8662) 


I ([Ca II] A7291) 


I ([Ca II] A7323) 


501 


15 


q 7 -l n 20 

i_> . 1 _1_ KJ . 


i An + 08 

-L.rtU _1_ KJ . UO 


n 40 -J- n 02 


n 41 + n (12 

KJ.'-tl. _l_ KJ • \J iJ 


505 


19 


q x. _|_ n 90 


i of. _|_ n 07 






508 


22 


3 2 4-015 

•>._ _1_ W.-LtJ 


i 1 s + 07 

-L.J-O _1_ W.L/l 






514 


28 


2.7 ± 0.15 








524 


38 


2.1 ± 0.12 


n sf; + n Ofi 

U.UU _1_ VJ.UiJ 


U.JI _1_ \J,\j±- 


n he + n ni 


530 


44 


1.6 ±0.12 








540 


54 


1.3 ±0.10 








554 


68 


1.2 ±0.10 








562 


76 






0.17 ±0.009 


0.15 ±0.009 


583 


97 


1.1 ±0.10 


0.45 ±0.05 


0.12 ±0.008 


0.10 ±0.009 


589 


103 


0.98 ± 0.09 








607 


121 


0.65 ± 0.07 


0.37 ±0.04 


0.051 ±0.003 


0.053 ± 0.004 


742 


256 


0.39 ± 0.05 


0.22 ±0.03 


0.022 ±0.001 


0.021 ±0.001 



a JD - 2454 000.00 
b Phase is in days with respect to the explosion date JD 2 454 486 ± 4. 

Table 11. Evolution of the three kinematic components (broad, v^; intermediate, vi; narrow, v n ) in Ha and Ca II (A8662) features. The velocity widths 
(FWHM) are in kms" 1 . 



JD a 


ph b 


v b (Ha) 


v; (Ha) 


v n (Ha) 


v b (Ca II A8662) 


Vi(Ca II A8662) 


v n (Ca II A8662) 


501 


15 


3110 ±300 


900 ± 90 


220 ± 30 








504 


19 


3020 ± 300 


875 ± 90 




3000 ± 350 


850 ± 90 


264 ± 30 


508 


22 


2910 ± 350 


870 ± 85 










514 


28 


2620 ± 370 


850 ± 90 










524 


38 








2280 ± 300 


753 ± 60 




530 


44 


2460 ± 250 


844 ± 90 


240 ± 30 








540 


54 


1940 ± 200 


790 ± 90 


260 ± 40 








554 


68 


1950 ± 250 


770 ± 95 


245 ± 30 








583 


97 


1920 ± 200 


750 ± 80 




2040 ± 200 


720 ± 80 




589 


103 


1610 ± 180 


680 ± 75 










607 


121 


1450 ± 100 


670 ± 90 










610 


124 








1340 ± 200 


660 ± 70 




669 


183 


1300 ± 90 


560 ± 70 




1290 ± 200 


550 ± 70 




742 


256 


1270 ± 100 


480 ± 70 




1260 ± 300 


540 ± 80 





a JD - 2454 000.00 
' Phase is in days with respect to the explosion date JD 2 454 486 ± 4. 



anism for line photons favoured by the s mall photospheric veloc- 
ity observed for SN 2008S at all phases jPessart et al1l2008l) . The 
electron scattering in a dense (Thomson optical depth of ~ 3- 
4) CS shell ly ing beyond the SN photosphere was suggested by 
Chugaij d200ir) as an explanation for the very broad wings of Ha 
in SN 1998S. We compare the Ha profiles of SN 1998S and SN 
2008S in Fig. [12] The redshift in the Ha peak decreases after about 
100 days (Fig. I16t . although the uncertainties in the wavelength 
calibration mean that this shift has only a modest significance. The 
EW(Hq) increases with time, becoming very large at late phases 
(about 900A). There are two processes proposed for the origin of 
H emission lines: photoionisation with subsequent emission by re- 
combination or collisional excitation jDrake & Ulrich|[l980l) . The 
first can not explain the increase of the EW(Hq) with time so colli- 
sional excitation of Ha from the n = 2 level seems to be the more 
likely mechanism to explain Balmer emission, especially at the lat- 
est phases. In the first spectra H/3 shows a multicomponent profile 
similar to Ha (Fig. |9), but this feature is too weak to analyse its 
profile in detail. The Ha/H/3 ratio is larger than the recombination 



value of about 3 and increased from 4 at 15 days to 10 at 256 days 
after the explosion. This ratio depends critically on several parame- 
ters: electron density (N e ), electron te mperature, external rad iation 
field strength and optical depth in th q (Dra ke & Ulrichll980h . The 
three most important processes which can change the line intensity 
with respect to the recombination value are: Balmer self absorption 
(important at low density), collisional excitation a nd de-excitation 
(important at high density) dDrake&Ulrlchll 19801) . 

Collisional excitation processes cause the Ha/H/3 intensity 
ratio to increase, and the N e value at which these processes become 
important inversely depends on th q - At the first epochs the Ha/H/3 
intensity ratio indicates of 10 10 cm -3 ^ N e ^ 10 12 cm -3 given 
the observed temperatures in SN 2008S (see Sect. |4) so the colli- 
sional processes may have an important role. The Balmer decre- 
ment, observed in SN 2008S, may be a sign of the high optical 
depth and the interaction with a high density CSM. In particular, its 
evolution is likely led by the increase of tho, with time. 



Calcium 

The strength of the [Ca II] doublet decreases by a factor of < 



10 
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Figure 9. Line profiles at about two weeks after the explosion epoch (JD 
2 454 486): in the top panel the multicomponent profile of Ha in the high 
resolution ISIS spectrum (the three Gaussian component are shown as thin 
dotted lines and their combination is shown as thick dotted line), in the mid- 
dle panel H/3 profile in the NOT spectrum, in the bottom panel the profiles 
of Fe II lines in the ISIS spectrum. Wavelength is in the observer frame. The 
flux densities have been normalised to the line peaks. 



Figure 11. Temporal evolution of Ha kinematic components and intensity. 
From the top to the bottom panel: evolution of the velocity of the narrow 
(v n ), intermediate (vi), broad (v^,) component, the integrated flux of Ha 
and the R band. Velocity widths (FWHM) are in kms —1 and integrated 
fluxes are in units of 10 -14 erg cm -2 s _1 . 



geometry in the CSM. 



over 100 days and line ratio Ha/[Ca II] ranges from 9 to 18 
over about 260 days (Table ITOt . We did not see any temporal 
evolution of the intensity ratio for the [Ca II] doublet (always 
around the value of 1±0.2) and of the position of the peak. The 
[Ca II] doublet very likely originates in a lo wer density region with 
respec t to the Ca IR triplet. As discussed bv lChevalier & Franssonl 
( 1994), the [Ca II] doublet originates from radiative de-excitation 
from the metastable 3d 2 D level, which is highly populated. The 
probabilities for the two [Ca II] transitions originating within the 
same multiplet are approximately equal and the threshold to ionize 
Ca + to Ca ++ from this metastable level is at 1218.8 A. The gap 
between this threshold and the Lya line is 3.2 A corresponding to 
about 800 kms . Only gas with higher velocity can ionize Ca 
II by absorbing Lya photons suppressing [Ca II] emission. The 
presence of strong [Ca II] doublet in the first SN 2008S spectra 
may be indicative of little Ca II being ionized, and hence a narrow 
Lya line. The [Ca II] doublet shows a velocity width (FWHM) 
of about 250±60 km s~ that is comparable to that of the narrow 
component of both Ha and Ca II triplet. This velocity width 
remains constant while the line profile shows an evolution: it is 
asymmetric during the first three months after the explosion and 
starts to be more symmetric after this epoch (Fig. 1 14b. In Sect. 17.21 
we propose an explanation for this behaviour based on a toroidal 



The Ca II NIR triplet must originate in a different region 
from that which produces the [Ca II] doublet since it shows both 
a multicomponent profile (Fig. 1151 ) and a decreasing redshift 
with time (Fig. I16t . Velocities for the different components 
are reported in Table QT] and are consistent with those of Ha. 
The intensity ratios of the infrared triplet to Ha and to [Ca II] 
increase with the time (Table UOt . To investigate if the peak 
position of Ha, the [Ca II] doublet and the Call NIR triplet 
show any temporal evolution, we again exploited the BIC factor. 
We compare the BIC values obtained by fitting the spectral 
time series with either a position constant in time or a temporal 
evolution (parameterized as a straight line with a non null slope), 
ABIC = BIC s io P c — BICconst • For Ha there is only a marginal 
evidence for an evolution, ABIC ~ 0.3. For the Ca II doublet, 
evolution is not favored (ABIC ~ —1.0 and —0.8 for the first and 
the second line, respectively). Finally, there is strong evidence for 
a decreasing redshift for the Ca II NIR triplet, ABIC ~ 13.9. 

Oxygen, Iron and Sodium 
The [O I] (AA6300,6364) doublet is visible in the first high res- 
olution spectrum with a very low velocity width (FWHM), about 
80 kms -1 . These are collisionally suppressed at high density 
so likely originate in a slow moving and low density region. O 
I(A8446) appears at the latest phases with a velocity width of about 
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Figure 12. Comparison between Ha profiles in SN 1998S (phase 9 days) 
and SN 2008S at two different phases (15 and 19 days). Phase is in days 
after the explosion epoch (JD 2 450 869 for SN 1998S and JD 2 454 486 
for SN 2008S). Wavelength is in the observer frame. A Gaussian profile that 
matches the line peak is shown as a dashed line. 



530 km s -1 similar to those of the Ha and Ca II intermediate com- 
ponents at the same phase. The Fe II lines, visible in the first high 
resolution spectrum, appear to have only a narrow component with 
velocity width around 200 kms _1 and the same profile as [Ca II]. 
The Na I D feature, visible in absorption in the early spectra, seems 
to show two blended lines: the Galactic contribution and the host 
galaxy doublet. In the spectra acquired at 182 and 256 days this 
feature clearly appears in emission. This evolution is indicative of 
the circumstellar origin of the Na I D and is a clear sign of the 
high density of CSM. The evolution of EW(Na I D), illustrated in 
Sect. 12.3 1 may due to an evolution of the ionization conditions in 
the CSM and in the ejecta of SN 2008S since the EW is related to 
the ionization stage of Na I. 



6 THE PROGENITOR OF SN2008S 

6.1 Optical and NIR pre explosion images 

We carried out an independent analysis on the deepest pre- 
explosion optical and NIR images t hat we could locate, similar to 
that presented in lPrieto et al.l d2008l) . We added some further data 
and recalculated all the limits as some values in that paper were 
taken from heterogeneous references rather than original data. We 
used optical images from the Gemini North and Large Binocular 
Telescope and NIR images from the Bok 2.3 m telescope. To en- 
sure an accurate positioning of SN 2008S on all the pre-explosion 
images we used an image of SN 2008S obtained with the Auxiliary 
Port Imager (AUX) on the WHT on 2008 February 08. The cam- 
era has a pixel scale of 0'.' 1 1/pixel, and a 1024 pixel TEK2 CCD. 
A single R band image of 300 s was taken, and although it was at 
high airmass (2.47) the image quality was l'.'03. The position of 
SN 2008S was determined on all the pre-expl osion images to an 
accuracy of 0'.' 1 and we confirm the results of IPrieto et al.l ( 2008) 
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Figure 13. The peak position and profile of Ho at at 19 days, about 2 
months, 4 months and 9 months after the explosion epoch (JD 2 454 486). 
Wavelength is in the observer frame. Flux density has been normalised to 
the line peak at each phase. 
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Figure 14. Profile of [Ca II] doublet at about one month, 3 months and 
4 months after the explosion epoch (JD 2 454 486). Wavelength is in the 
observer frame. Flux density has been normalised to the line peak at each 
phase. 
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Figure 15. Profile of Ca II NIR triplet at about 20 days, one month, 3 
months, 4 months and 9 months after the explosion epoch (JD 2 454 486). 
Wavelength is in the observer frame. Flux density has been normalised to 
the line peak at each phase. 
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Figure 16. The peak position of Ha, [Ca II] (A7324) and Ca II (A8662) as a 
function of phase. Phase is in days after the explosion epoch (JD 2 454 486) 
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that there is no optical or NIR counterpart detected. We calculated 
the 3a detection limits for each of the optical images from the LBT 

and Gemini images and these are reported in Table 1 12] 

The K' band image was obtained from iKnapen et all (2003) 
with the Bok 2.3 m. telescope of the Steward Observatory on 1999 
October 17 with the PISCES camera, a HAWAII HgCdTe array of 
1024x 1024 pixels of 0'.'5 on the sky. The dithered exposures re- 
sulted in co-added total exposure time of between 1000-2000 sec- 
onds, depending on field location, and the image quality was 2 / .'3. 
The WHT AUX R band image was used to locate the position of 
SN 2008S on the frame by matching the positions of 12 stars lo- 
cated within 1.7' of the SN in both R and K images. The geomet- 
ric transformation resulted in an RMS to the fit of 108 milli-arcsec 
and at this transformed position there is no detection of a source 
in the K' band image. The zero-point for this frame was estimated 
using five 2MASS stars located close to the SN and the calculated 
3a limiting magnitude was estimated to be K = 18.0 mag. The 
mean magnitudes of several of the faintest point sources visible in 
the vicinity of the SN was K — 18.3 ± 0.3 mag. Hence we adopt 
if = 18 mag as the sensitivity limit of this frame. 



6.2 Spitzer MIR pre explosion images 

As discussed in IPrieto et"al I d2008l) . several epochs of archival 
MIR imaging from the Spitzer Space Telescope are available. We 
analaysed IRAC (3.6-8.0 /im) images with the longest integration 
times from programmes 3249, 20256, and 30292 (P.I. Meikle). 
Aperture photometry was carried out on the post-BCD images us- 
ing GAIA. A circular aperture of 2'.'0 was used. Aperture correc- 
tions were derived using the point response function frames avail- 
able from the Spitzer Science Center. The residual background 
level was measured using a clipped-mean sky estimator, and a con- 
centric sky annulus having inner and outer radii of 1.5 and 2 times 
the aperture radius, respectively. The resulting flux densities are 
listed in Ta ble [T2l and are consistent with the values reported in 
IPrieto et al.1 J2008h . We used the AUX image of SN 2008S and the 
deep, wide field Gemini i band image to determine the position of 
SN 2008S on the Spitzer 4.5, 5.8 and 8.0 pirn images using differen- 
tial astrometry. In each case the uncertainty in the positioning of SN 
2008S on the Spitzer images was / .'3. This uncertainty is a combi- 
nation in quadrature of the geometric transformation RMS and the 
uncertainty in the measurement of th e centroid of the progenitor 
object detected by IPrieto et al. (2008). The difference between SN 
2008S and the progenitor object is 07 16±0730 and 07 14±0730 for 
the 5.8 and 8.0 ^im images respectively. However, there is a small 
difference in the position of the progenitor source in the 4.5 fim im- 
age. This is separated from the SN 2008S position by 0'.'51±0'.'33. 
On close inspection of the images it appears that the source in 
the 4.5 /im image may well be extende d compared t o its 5.6 and 
8.0 counterparts (see also Fig. 1 in IPrieto et alj «|2008)). This 
might suggest that the 4.5 fim source is actually a blend of two or 
more sources. It seems fairly secure that SN 2008S is coincident 
with the 5.8 and 8.0 sources and we have no evidence to sug- 
gest that they are extended or non-stellar. The suggestion that the 
4.5 [im source may be slightly extended could just be due to con- 
tamination from nearby, but unrelated flux at this wavelength. But 
the possibility remains that the source is a blend and this needs to 
be clarified with late, deep imaging with Spitzer and ground-based 
high r esolution s t udies. 

IPrieto et alj (2008) also reported aperture photometry for the 
individual pre-explosion Spitzer epochs. They found no evidence 
for variability in the flux of the pre-explosion source over the 
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Figure 17. Pre-explosion images of SN 2008S. (a) The Gemini GMOS i' band image showing no progenitor detection to deep limits, (b) K band image from 
the Bok telescope, again showing no detection in the NIR. (c) & (d) : Spitzer images at 4.5 and 5.8 fim, originally presented by Prieto et al (2008). We find 
some evidence that the 4.5 fim may be slightly extended. This can only be tested with higher spatial resolution images once the SN fades. 



~ 1000 days covered. Our photometry is based on the combined 
Spitzer images. To detect possible variability on the source flux 
we compared the Post-BCD 4.5 fim images from 2004-11-25 and 
2006-08-12 where the pre-explosion source was well detected. We 
first aligned the images using the centroid coordinates of 13 point- 
like sources around the pre-explosion source position. The a ligned 
image s were matched and subtracted using the ISI S 2.2 dAlard 
2000) image subtraction package (for details see also lMeikle et al. 
(2006). No residual above the noise was apparent in the subtracted 
image at the pre-explosion source position. This confirms the find- 



ings of IPrieto et alj d2008l) on the lack of variability of the pre- 
explosion source. 



6.3 Analysis of the pre-explosion images 

To determine whether the precursor object detected in the Spitzer 
images is a viable stellar progenitor we modelled the MIR emission 
using the radiative transfer model DUSTY, in a manner similar to 
that appl ied to dusty red superg iants (RSG) and AGB stars in the 
LMC by Ivan Loon et al.1 ( 120051) . In Fig. [18] we show two fits. In 
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Table 12. Fluxes and magnitudes measured in the position of SN 2008S on 
pre-explosion images. 



Waveband 


Instrument 


Flux or magnitude 


Units 


u 


LBT 


> 25.1 


Vega Mag 


B 


LBT 


> 24.5 


Vega Mag 


V 


LBT 


> 24.5 


Vega Mag 


i 


GMOS-N 


> 24.4 


Vega mag 


K' 


PISCES 


> 18 


Vega Mag 


3.6/tm 


IRAC 


< 3.6 




4.5/tm 


IRAC 


21.3 ± 1.5 


idy 


5.8/im 


IRAC 


45.6 ±2.4 




8.0/im 


IRAC 


59.5 ±4.3 


/iJy 



warm dust 
interstellar 



the first case (model A) we show a warm circumstellar, spherical, 
dust shell with T dus t = 800 K (at the inner boundary), which has 
an optical depth of tv = 150 (rs^m ~ 2). We assume that the 
central exciting source is a blackbody of 3000 K (a cooler temper- 
ature does not make a major difference to the MIR output). The 
outer radius of the shell is 454 AU (2.2 milli-pc), but its thickness 
needs to be R oute r/Ri nner < 5 to fit the SED. The dust in this case 
consists of silicates with a typical grain size distribution. The in- 
ferred luminosity of the central star is logL/L© ~ 4.6 and with 
such a dense shell all of the stellar energy is absorbed by the dust 
and hence all appears as reprocessed MIR flux. While this may 
seem a plausible model and scenario for a dust shell around a mas- 
sive, embedded star, the extinction toward the central star would 
be Av = 163 mag. Hence in this case of a spherical, dense (but 
radially thin) circumstellar shell the optic al and K' bands off er no 
meaningful constraints on the stellar SED. IPrieto et alj d2008l) sug- 
gested a thermally radiating sphere to account for the observed SED 
(with a blackbody temperature Td us t = 440 K, logL/L© ~ 4.5, 
andR BB = 150 AU). 

The second case (model B) shown in Fig. [T8] is a block of 
optically thick interstellar dust with Av = 140 mag. The extinction 
law from lCardelli etalj dl989l) is used, although this does not fit 
particularly well around the 8 /im detection. The extremely high 
value of extinction would imply densities similar to those in the 
cores of molecular clouds and the central stellar source would need 
to be logL/L© ~ 6.9 (again we have assumed a blackbody of 
3000 K). 

The main problem with either of the above two scenarios is 
that the extinction we see toward SN 2008S can be robustly esti- 
mated to be less than Ay < 4 mag and is more likely to be of 
order Av = 2.2 mag. Hence if a spherical dust shell surround- 
ing an exciting source is the explanation for the precursor SED, 
then the explosio n must have destroyed nearly all the dust within 
the 45 4 AU shell. iBode & Evanl dl98Cl) . IWrightl dl980h and lDwekl 
dl983l) have suggested that the UV-optical luminosity from SNe 
could destroy dust within a spherical cavity around the progeni- 
tor star. Additio nally UV flashes h a ve been recently observ ed from 
type II-P SNe dGezari et ai] |2008l : ISchawinski etai] l2008). How- 
ever, these calculations assume a small optical depth, and a more 
detailed estimate of an optically thick CSM which is evaporated by 
UV-optical photons within a gamma-r ay burst beam has been un- 
dertaken by IWaxman & Draind (2000). They suggest that a GRB 
could clear very high column densities of dust (e.g. Av ~ 30 mag) 
within several parsec of the explosion. Taking these calculations, 
scaling the destruction radius R d by y/h ( where for a SN X-r ay/UV 
flash we take E^ 10 46 erg, as found in iGezari et al. I d2008l) and 
ISchaw inski et al] d2008h ) and reducing the effective flux by the 
solid angle of the GRB beam in comparison to the isotropic SN 




wavelength [^m] 

Figure 18. The MIR SED of the SN 2008S progenitor star and the two fits 
obtained using the radiative transfer model DUSTY. Model A (thin line) 
consists of a warm circumstellar, spherical, dust shell with T^nst = 800 K 
(at the inner boundary), which has an optical depth of Ty = 150. Model 
B (dotted line) consists of a block of optically thick interstellar dust with 
Ay = 140 mag. 



energy, a destruction ra dius of a few milli-parsec is found. In the 
optically thick regime of lWaxman & D rain el d2OO0l) the destruction 
radius does not depend strongly on the density, or Av, with Rd 
changing by a factor of less than 2 when nH changes by a factor 
10°. Hence it may be feasible for a dust embedded SN to clear a 
large enough cavity to become relatively unobscured, even at rela- 
tivel y high dust densitie s . 

IPrieto et"al] d2008h proposed that the properties of the precur- 
sor source are very similar to those of dust enshrouded R SG and 
AGB stars in the L MC as studied by Ivan Loon et all d2005l) . All of 
the carbon stars in Ivan Loon et alj d2005l) hav e log L/L© < 4.09 
which is somewhat too low to result in a CCSN. ISmarttetal.l f2009) 
show the lowest luminosity progenitors to be around 4.3 dex and 
stellar evolutionary models would suggest progenitors of luminos- 
ity < 4.1 dex would have initial masses of 5M0or less. Hence 
these C type stars are not plausible counterparts to the SN 2008S 
progenitor and are not viable progenitors for a CCSN. More lumi- 
nous and massive stars than these AGB objects (which are above 
the core-collapse threshold) tend to have less optically thick en- 
velopes. For example, IRAS 04516-6902 is likely to have Av ~ 
13 mag but would be too faint to match the precursor of SN 2008S. 
In fact all of the RSGs in LMC with luminosities above 4.3 dex 
would not be bright enough at 8 /im to account for the observed 
flux before the explosion of SN 2008S. In addition, many of them 
would be too bright in the K ox I band to be consistent with our 
observed upper limits (although adjusting the extinction could help 
hide this shorter wavelen gth flux). 



IRAS 05280-6910 dvan Loon etal]|2005h is an extreme case 
where the star is hardly detectable in the optical and might ap- 
pear an appealing source for comparison. However in this case the 
24 /xm flux of IRAS 05280-6910 would be much brighter than the 
limit set for the pre-explosion source. Hence we agree that the sug- 
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Figure 19. An HR diagram showing the positions of RSGs in the Galaxy 
from Levesque et al. (2005) (points) and the region (shaded region) in which 
red supergiant progenitors of normal II-P SNe have been seen (Smartt et 
al. 2008). The luminosity of the progenitor of SN 2008S is estimated as 
logL/L = 4.6 ± 0.3, and T cff = 3000 K is consistent with the SED 
(black point with error bar). However the T c ff is unconstrained from the 
DUSTY model. The dotted lines limit the region of luminosities allowed 
for the progenitor. The tracks are the STARS models from Eldridge & Tout 
(2004). 



sistent with the position we would expect in the HRD for a 6-8 Mq, 
star which has gone through 2nd dredge-up and has developed an 
O-Ne-Mg core within which electron-capture collapse could oc- 
cur. The luminosity of the progeni tor has been interpreted as an 
indication of a mass of 10-20 Mm jPrieto et alJbOQSl : ISmith et al.l 
I2008L iBond et al.ll2009l : feerger et alj|2009h . However, as seen in 
Fig. [19] the luminosities of 6-8 Mq progenitors after 2nd-dredge up 
are consistent with the bolometric luminosity of the progenitor of 
SN 2008S. 

So far our analysis has assumed that the pre-explosion MIR 
source is a single object and not, for example, an embedded cluster 
of stars within which the progenitor arose. One cannot definitively 
rule out the latter and we noted above that there is some evidence 
to suggest that the 4.5 /im image is extended. However, we would 
argue that it is unlikely to be a cluster for two reasons. Firstly, the 
total luminosity is not unusually high for a stellar source, and if it 
is a cluster that hosted a SN explosion then the progenitor is likely 
the dominant source of flux. Secondly, the cluster would still need 
to have an extinction similar to that found for the stellar source 
assumption and the dust would, presumably, be extended over a 
few parsec (typical cluster size). However, we estimated the likely 
destruction radius at only a few milli-parsec which seems too small 
to have an embedded object within a cluster clear the line of sight 
of intervening dust. The CS dust shell appears to be more plausible. 



gestion of lPrieto et al. I d2008l) of a dust enshrouded RSG is initially 
very appealing, and the dust destruction even at such high column 
densities does not appear unrealistic. However, we cannot easily 
match SN 2008S progenitor quantitatively with any of the known 
dusty RSG or AGB stars in the LMC and Galactic samples, in full 
agreement across the optical and NI R non-detections. Our c onclu- 
sions are consistent with the work of [Thompson et all (2008). who 
show that stars with the same MIR properties as the progenitor of 
SN 2008S are very rare in the nearby spiral M 33. They find only 
~ 10 objects with similar magnitudes and colour which they asso- 
ciate with the high luminosity tail of the AGB sequence. They con- 
clude that this phase is a short period in the lives of a reasonable 
fraction of massive stars, perhaps linked to large mass ejections and 
subsequent dust formation episodes in the last ~0.1% of a stellar 
lifetime. The fact that these types of stars are not common in the 
LMC studied populations supports this conclusion. 

In Fig. [19] we show an Hertzsprung Russell diagra m (HRD) 
with the STARS model tracks of lEldridge & Toutl d2004h. The po- 
sitions of Galactic RSGs are shown (from Levesque et al.l 120051) 
along with the region in which RSG pro genitors of recent n earby 
type II-P SNe lie (logL/L Q = 4.3+g j ISmartt et alj|2009h . The 
luminosity of the progenitor has been estimated at logL/L© = 
4.6 ± 0.3. Although we have estimated the temperature of the ex- 
citing source as T e g = 3000 K in the DUSTY model, this num- 
ber is not well constrained, as a hotter star can be placed inside a 
denser, geometrically thinner envelope and produce a similar SED. 
The dotted lines at the extrema of the luminosity ranges show where 
the progenitor may lie. The tracks of 6-8 Mq stars which have gone 
through 2nd d r edge-u p are shown as dotted line. As discussed by 
Eldrid ge"et alj |2007), many different stellar evolutionary models 
predict that super- AGB stars (in the 6-9 Mq range) that have gone 
through 2nd dredge-up can rise to higher luminosities and lower 
effective temperatures than their higher mass counterparts. These 
stars may be thermally unstable, pulsating, and prone to large mass 
ejection events. The luminosity of the SN 2008S progenitor star is, 
within the uncertainties of both the measurements and models, con- 



7 THE NATURE OF SN 2008S 
7.1 Pre-explosion CSM 

Super-AGB stars lose mass at a high rate and their circumstellar 
medium have a high density and complex geometry. The circum- 
stellar environment of SN 2008S suggests such a progenitor star. 
The progenitor analysis demonstrates the presence of a dusty opti- 
cally thick shell around SN 2008S with an inner radius of nearly 
90 AU , while the MIR echo analysis indicates a second, dusty 
outer shell of inner radius ~2000AU (Fig. 120b. The two different 
shells imply that the progenitor star experienced a variation in the 
mass-loss rate. The physical size of the inner shell indicates that 
it was probably circumstellar and the result of m ass loss due to a 
steady wind rather than eruptive ejection (see also Thomps on et al.l 
(2008)). If we assume that the optical depth to the source scales as 
r~ 1 in a freely expanding wind, tv to the progenitor was about a 
factor of ~ 10-20 larger before explosion i.e. the inner shell was 
highly-obscuring. The dust within this shell was probably evapo- 
rated in the explosion. In contrast the much larger distance of the 
dust in the outer shell meant that it was not destroyed by the ra- 
diation from the SN photosphere. Instead, the enormous MIR flux 
observed at 17.3 days can be explained as an IR echo from the outer 
shell dust whose mass we estimate to be ~ 10~ 3 Mq. We also pro- 
pose that, if we invoke a toroidal geometry in the outer shell, then 
the asymmetric profile of the [Call] and Fe lines can also be ex- 
plained. This is addressed below. 



7.2 The photometric and spectroscopic evolution 

SN 2008S exhibited a very slow photometric evolution and al- 
most no spectral variability during the first nine mo nths, implying a 
long photo-diffusion time and a high density CSM dSchlegelll 990). 
Nevertheless, photometric and spectroscopic observations suggest 
that there are three distinct phases: the maximum light phase (0- 
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50 days after the explosion), the flattening of Ha intensity (60-100 
days) and the NIR excess phase (from about 120 days). 

During the first two phases the SN 2008S optical-NIR SED 
can be fitted with a single blackbody with a radius of about 
2 x 10 14 cm. The blackbody temperature and radius declined mono- 
tonically, with the temperature falling from ~ 8000 K to 5000 K. 

During the first phase, all the spectra showed prominent emis- 
sion lines of Ha, the [Ca II] doublet, the Ca II NIR triplet, faint Fe 
II emission lines and Na I D in absorption. None of the H, Ca or Fe 
lines exhibit P-Cygni profiles, providing additional evidence of a 
high density CSM. Some of the Na I D absorption probably has an 
origin in the CSM, given the observed high density and rapid evo- 
lution of the EW during the first month. The only forbidden lines 
visible in the early spectra are the intense [Ca II] and weak [O I] 
doublets, probably produced in the low-density inter-shell zone. 

At about 60 days after the explosion the decline of SN 2008S 
light curves steepened, while the Ha intensity remained constant 
until about 100 days. During this phase the Ha narrow component 
disappeared and the [Ca II] profile lost its earlier asymmetry. We 
suggest, as an explanation of the change in the [Ca II] profile, that 
the progenitor CSM took the form of a toroid, viewed almost edge- 
on. At early times, obscuration of the [Ca II] emission zone by the 
outer shell dust would affect predominantly the fraction of the gas 
moving towards us, thus attenuating the blue wing. However, by 
60 days, much of the ejecta would have reached the inner shell. (It 
would take a velocity of only ~ 2500 kms -1 .) This would be ex- 
pected to produce a strong burst of hard radiation arising from the 
ejecta-shell impact. This, in turn, might evaporate the outer shell 
dust resulting in the increasing symmetry of the [Ca II] line pro- 
files. However, there are several difficulties with this scenario, (a) 
We surmise that the [O I] doublet formed in the same low density 
region as did the [Ca II] feature, and yet no comparable wavelength 
shift is seen in the [O I] feature, (b) The IR echo model indicates an 
extinction through the outer shell dust of only about A = 0.15 mag 
at 7300 A. Even allowing for the toroidal geometry, it is not clear 
that this would be sufficient to account for the profile change, (c) 
If part of the ejecta/inner CSM impact radiation flowed inward, it 
might destroy or inhibit the formation of new ejecta dust. However, 
it may be that the impact radiation would be severely attenuated 
by intervening material before it could reach the dust-formation 
zone. Further discussion of the [Ca II] profile evolution is beyond 
the scope of this paper. We note that the release in 2009 August of 
the Spitzer observations of SN 2008S at ~ 180 days will provide 
a useful test of the outer shell dust evaporation scenario. We con- 
clude that while the NIR excess provides clear evidence for newly- 
formed dust, the evolution of the [Call] feature is somewhat more 
difficult to explain. 

The last phase started after about 120 days when the appear- 
ance of the NIR excess indicated the presence of an additional, 
warm component. The temperature and radius (and therefore the 
luminosity) of the hot component showed a slower decline dur- 
ing this phase. The warm component cooled from 1400 K at 160 
days to 1200 K at about 300 days, while its radius and luminos- 
ity increased. An explanation for the NIR excess is thermal emis- 
sion from newly-formed dust in the ejecta or in a cool dense shell 
formed by the ejecta-inner CSM impact. At this epoch the decline 
rate of the SN 2008S light curves flattened to 1.3 to mag/lOOd, de- 
pending on the band. The bolometric light curve shows a decay rate 
of 0.88 mag/lOOd, very similar to that of 56 Co. The velocity width 
of the broad component and the intensity of Ha declined slower 
while the velocity width of the intermediate component decreased 
as in the first phases. The [Ca II]/Ha intensity ratio halved by 260 
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Figure 20. Schematic illustration of the pre-explosion geometry of 
SN2008S. 



days while the (Ca II triplet)/Ha declined only slightly. During this 
later phase the Na I D and O I appeared in emission. 



7.3 Comparison with NGC 300 OT2008-1 and M 85 
OT2006-1 

Several recent papers have suggested an analog y between SN 
2008S and two other tr a nsients : M 85 OT2006-1 dKulkarni et ail 
l2007t IPastorello et ail l2007h and NGC 300 OT2008-1 



(Tho mpson et al.l 1200 



iBerger et alj 12009] ; iBond et alj l2009h . 
The progenitors of the SN 2008S and NGC 300 OT2008-1 
transient s have been dete c ted only on Spitzer pre-explosion 
image s jPrieto et al] l2008t iThompson et alj 120081 ; IBerger et all 
120091) which indicate s very similar p r operti es and geometry of 
the obscuring CSM. Thompson et alj d2008h has proposed that 
these transients are members of a new class given the similar 
pre-explosion properties (see Sect. [63). The extensive multi- 



wavelength monito ring campaigns for both these transients ( 
ISmithet alJ d2008l) and this p aper for SN 2008S; IBond et al.1 
(2009) and IBerger etafl d2009l) for NGC 300 OT2008-1) show 
that they are indeed very similar in their kinetic and radiative 
energies. M 85 OT2006-1 had a peak R band absolute magnitude 
of ~ — 12 , a peak luminosity L p ~ 2 x 10 40 ergs _1 and a 
total radiated energy in the first two months E p h ~ 6 x 10 46 erg 
dKulkarni et all 120071) , while NGC 300 OT2008-1 had a V band 
peak of —13, a luminosity of 1.6 x 10 40 erg s _1 and a total energy 
of ~ 2x 10 47 erg and no radio or X -ray emission as SN 2008S 
dBond et al]2009l ; IBerger et alj|2009h . 

The light curves of both transients are very similar to that of 
SN 200 8S although in the la test phases M 85 OT2006-1 faded 
quickly 1 Kulkar ni et alj 120071) . as shown in Fig. [2J] In the first 
phases of its evolution M 85 OT2006-1 showed a warm infrared- 
bright component likely due to an IR echo by circumstellar dust 
surviving the explosion dPrieto etalj|2008h . The spectra of both 
M 85 OT2006-1 and NGC 300 OT2008-1 s howed B aimer lines 
the [Ca II] doublet and the Ca II NIR triplet dKulkarni et aHl2007l; 
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IPastorello et alj|2007l : lBond et alj|2009l : IBerger et aI1l2009ri . More- 
over, in NGC 300 OT2008-1 He I (in emission), and Ca H&K and 
OI A8446 (in absorption) are visible while M 85 OT2006-1 showed 
also prominent K I lines. The velocity width of Ha is very similar 
in SN 2008S and NGC 300 OT2008-1, while in M 85 OT2006-1 is 
narr ower. 

ISmith et al] fc008b have interpreted SN 2008S as a SN impos- 
tor analogous to the eruptions of LBVs. They proposed that SN 
2008S was a super Eddington outburst of a star of about 20 Mq, 
highly obscured because an outburst which had occurred shortly 
after the recent blue loop transition. However, a 20 M© star has a lu- 
minosity of log L/L0 ~ 5.0 — 5.3 which is not consistent with the 
total MIR luminosity of the progenitor of log L /Lq ~ 4.6 derived 
from the analysis of the pre-explosion images. [Smith et all (2008) 
also suggested that the spectral similarity of SN 2008S in outburst 
with the hypergiant IRC+10420 might indicate that SN 2008S was 
also a star of similar evolutionary state. However, the similarity be- 
tween the SN 2008S and IRC+10420 spectra only points to similar 
physical conditions of the regions where the emission lines form. 
The fact that the T c ff of the continuum of SN 2008S decreases dra- 
matically, while the emission-line spectrum does not evolve, sug- 
gests that the lines are not formed in a stellar like expanding pho- 
tosphere. IRC+10420, in its quiescent phase, has a luminosity 20 
times higher than that of the MIR progenitor but 50 times lower 
than that of SN 2008S at maximum. Finally IRC+10420 is not en- 
shrouded by a dusty shell and likely has a different mass loss his- 
tory with respect to that of progenitors of these transients. Overall, a 
physical or evolutionary link between possible LBV-like outbursts 
and SN 2008S do es not seem convincin g to us. 

In contrast iThompson et al] ( 120081) compared the MIR prop- 
erties of known LBVs in M 33 to the progenitors of SN 2008S 
and NGC 300 OT2008-1 and claimed that a LBV explanation is 
unlikely for these transients since the LBV luminosity is higher 
and the LBV MIR colours are much bluer than those of the two 
transients. Moreover, they stressed that the time scale of the LBV 
variability is not consistent with the lack of variability of the MIR 
progenitors. We would agree that this is an evidence against the 
stellar eruption scenario, at least in any LBV or LBV-l ike event. 

As an explanation for NGC 300 OT2008-1, iBerger et al] 
d2009h and lBond et aT] d2009r) also favour a stel lar eruption which 
is not unlike that proposed bv lSmith et al] {2008). But nevertheless, 
the issues discussed above still argue against this interpretation for 
SN 2008S. 

iBond et al] J2009b proposed as the SN 2008S-like transient 
progenitors heavily dust enshrouded luminous stars of about 10- 
15 Mq, likely an OH/IR sources, which have begun to evolve on 
a blue loop toward warmer temperatures. During this transition 
progenitor stars reached a state in which they exceeded the Ed- 
dington limit for their luminosities and masses and suddenly initi- 
ated outflo ws. The reason of t hese eruptions nevertheless remains 
uncertain. IBerger et al] j2009h explained the nature of NGC 300 
OT2008-1 and SN 2008S as an eruption of a blue supergiant or 
pre-WR star with a mass of 10-20 Mq likely in a binary system. 

7.4 Comparison with SN 1998S and SN 1979C 

The optical light curves of SN 2008S are surprisingly similar to 
those of type Iln SN 1998S and type ILL SN 1979C (Fig. |22j, 
although these we re much more luminous (Mb = —19 . 6 from 
iFassia et al] feOOCb and M fl = -19.4 from lPanagia et al] dl98Ch . 
respectively). In the NIR bands the decline rates of SN 2008S at 
~ 170 days after the explosion are similar to those of slowly de- 
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Figure 21. Comparison between light curves of the transients SN 2008S, 
M85 OT2006-1 iKulkarni et all |2007|) . NG C 300 OT2008-1 iBond et al] 
120091) , t he subluminous type II-P SN 1999br iPastorello et alj2004) and SN 
2005cs IPastorello et all200dlTsvetkov et alj200dlPastorello et all2009» . 
and the SN impostor SN 1997bs IVan Dvk et alj2000h~ 

dining CC SNe templates of iMattila & Meikld d200lh based on 
SN 1979C and SN 1998S at the same epoch (J ~ 0.9mag/100d, 
H ~ 0.3mag/100d and K ~ Omag/lOOd). In Fig.|22]the R band 
and UBVRIJHK quasi-bolometric light curves of SN 2008S are 
compared with those of SN 1998S, while in Fig.|23]the evolution of 
the B - V, V - R and V - K colours is shown for SN 2008S and 
SN 1998S. The result of this comparison is intriguing: the overall 
photometric evolution of these events is very similar, the only two 
differences being the absolute luminosity and a broader peak for 
SN 2008S. 

Hydrodynamical modelling of type ILL SNe has suggested 
three quite diffe rent evoluti onary scenarios and explosion mech- 
anisms. ISwartz et al] ( 1 199 ll) modeled the collapse of a O-Ne-Mg 
star, suggesting the electron capture as a mechanism for ILL SNe. 
lBlinnikov & Bartunovj ill 9931) proposed a large supergiant progeni- 
tor. while lYoung et al] d2005h presented a two component model of 
a GRB afterglow with underlying SN ejecta. However, there is still 
not a clear consensus on the bright type Iln and II-L events such as 
SNe 1998S and 1979C. 

The spectral evolution of SN 1979C and SN 1998S is quite 
different with respect to that of SN 2008S but there are some 
common characteristics. Both SN 1979C and SN 1998S showed 
strong emission lines and absenc e of P-Cygni profi l es during the 
first m onths after the discovery dFassia et al] |2000| : iBranch et al] 
1 198 lb . The high density of the CS shell in the case of SN 1998S 
implies that a cool dense shell (CDS) forms at the interface of SN 
ejecta and the CSM. The absence of bro ad P-Cygni pr o files was ex- 
plaine d with obscuration by the CDS (Chugai 2001). Fass ia et all 
(2001) showed that absorption troughs appeared in the spectra of 
SN 1998S after only 12 days, and suggested that the SN ejecta had 
over-run the inner CS by that point. 

SN 1998S and SN 1979C have been suggested to be the results 
of explosions of a red supergiant with an extended envelope, R~ 
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Figure 22. Comparison betwe en light curves of SN 2008S (filled s ymbols). 
SN 1998S (empty trian gles) <Liu et alj|200(i [Fassia et afll200ll) and SN 
1979C (empty squ ares) (Balinskaia et al. 1980; de Vaucouleu rs et aljl98ll ; 
iBarbon et all 19821) : in the top panel the R band light curves (the SN 1998S 
and SN 1979C light curves are both shifted by 4 mag) and tin the bot- 
tom panel the quasi bolometric light curves of SN 2008S and SN 1998S. 
Phase is in days after the explosion epoch (JD 2 450 869 for SN 1998S, JD 
2 443 979 for SN 1979C and JD 2 454 486 for SN 2008S). 
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Figure 23. Comparison between V — K, V — R and B — V colours of SN 
2008S (filled symbols) and SN 1998S (empty symbols). The colours of SN 
1998S are corrected for Ag Q a i = 0.86 and Ag^ nt = 0.09. Phase is in 
days after the explosion epoch (JD 2 450 869 for SN 1998S, JD 2 443 979 
for SN 1979C and JD 2 454 486 for SN 2008S). 



(l-10) x 10 3 Rq, and a moderate mass of ejecta ~ 5M Q .[chuHai 
suggested also that the origin of the CS shell around SN 
1998S may be a violent mass loss during the Ne and O burning in 
cores of ~ 1 1 M e . The mass-loss rate of the progenitor by optical, 
radio and X ray estimates, is about 1 — 5 x 10 -4 Mq yr _1 . Likely 
the wind of the SN 1998S and SN 1979C progenitor underwent 
significant changes in the rate of mass outflow resulting in several 
distinct CSM shells. 



7.5 SN 2008S: an electron capture SN 

An alternative explanation of SN 2008S is the explosion of a mas- 
sive star in a low, or moderate, energetic SN. We show that the 
observed characteristics of SN 2008S and its progenitor share char- 
acteristics that are common to models of electron-capture super- 
novae (ECSNe). Such models have been developed extensively in 
the last three decades. Thomps on et alj d2008h suggested that the 
ECSN was a viable explanation for these SN 2008S-like objects 
based on the progenitor properties. The energetics of the events and 
their evolution now add weight to this argument. The strongest ev- 
idence for this interpretation is based on the observational discov- 
ery that the light curve tail of SN 2008S follows the decay rate of 
56 Co. The observation of 56 Ni and 56 Co decay can be used to dis- 
tinguish between a SN explosion and multiple shell interaction sce- 



narios jDessart et alj20 08). The tail phase luminosity decay rate is 
a strong argument that 56 Ni was produced by SN 2008S, and there 
is no physical channel to produce this radioactive isotope other than 
explosive burning of oxygen and silicon at temperatures in excess 
of ~ 10 9 K. 

A star in the mass range of ~ 8-10 Mq can form an electron 
degenerate O-Ne -Mg core during the final stage , as it becomes a 
super - AGB star dRitossa et ai]|l999l ; Isiess 2007; lEldridge & Tout! 
2004). Super-AGB stars may end their lives either a s massive O- 
Ne-Mg white dwarfs or as ECSN before Ne ignition dMivaii et al.l 
198CA Nomot cl 1 19841; Iffillebrandt et al] 1 19841; iBaron et alj 1 1987 



Mavle & Wiisonll 19881 ; IPoelarends etai]|2008l ; IWanaio et al]|2003l 
20081) . The final fate of super-AGB stars depends on the com- 



peting effects of core growth and mass loss during the late evo- 
lutionary stages, particularly during carbon burning. If the core 
mass reaches the Chandrasekhar mass, high pressure and density 
lead to electron capture onto 24 Mg and 20 Ne, the electron de- 



iron core is formed l 


Mivaii etai]|l98fj(: iMiyaii & Nomotoll 19871; 


Hashimoto etal J 1993 


;|Poelarends et al. 2008; Wanaio et al.l2008). 



If the mass-loss rate is high enough the envelope is lost before the 
core reaches the Chandrasekhar mass and the star ends its life as an 
O-N e-Mg wh i te dwarf. 

lNomotol ( ll984ll987t) estimated that stars with a mass of 8- 
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10 Mq can become ECSNe but recent work has su ggested a more 
narrow mass range ( Siess 2007; Poelarends et al. 2008) and a lower 
mass limit of ~ 9 Mq . However. IPodsiadlowski et al] d2004l) sug- 
gested that the initial mass range may be wider if one considers 
binary systems, about 8-11 Mq. ECSNe and their progenitors are 
predicted to show three properties that might allow us to distin- 
guish them from ordinary CCSNe: they might produce low energy 
explosions; the enormous mass-loss rate of the progenitor star in 
the super-AGB phase may produce signatures of CSM interaction 
in SN the light curve; the ECSN progenitors have l uminosities of 
the or der of 10 5 Lq and cool effective temperatures dEldridge et al.l 
l2007h . All these characteristics appear consistent with SN 2008S 
and its progenitor star. An important implication of the low ex- 
plosion energy is a small 56 Ni mass produced by the explosion. 
In ECSN explosion models, iKitaura et alj |2006) estimated an ex- 
plosion energy of 1-2 x 10 50 erg, a Ni mass of about 10 -2 Mq, 
and an ejecta velocity after the shock breakout from the stellar 
surface of ^ 3000 kms -1 . The estimated values of energy and 
56 Ni mass in SN 2008S are somewhat lower than the prediction 
from Kitaura et alj d2006l) . but more recent models of ECSNe by 



Wana io et al 



d2008h the observed properties of SN 

2008S. These weak explosions give an ejected 5B Ni mass of 0.002- 
0.004 M . 

The progenitor of SN 2008S is too faint to be a luminous RSG, 
and too bright to be an intermediate-mass AGB star. It could thus 
be a super-AGB star or a "low-mass" (close to ~ 8 Mo also) RSG. 
One argument in favour of a super-AGB star progenitor, is that low- 
mass RSGs generally do not lose mass fast enough to avoid explod- 
ing as a RSG, by far. So they still have massive mantles, but fairly 
diluted circumstellar envelopes. AGB stars and super-AGB stars 
lose mass at high rate for longer prior to their ends, which in case 
of AGB stars means they lose their mantles before core-collapse 
could occur whereas in super-AGB stars it's a race between mass 
loss and core growth that determines who wins. The stellar mantle 
prior to explosion would be thin anyway, one would expect, and 
the circumstellar medium would be dense. This strong mass loss is 
an obvious way to account for th e optically thick shel l arou nd SN 
2008S and NGC 300 OT2008-1. Weaver & Woosievl Jl979h sug- 
gested that the super-AGB stars might produce strong flashes in the 
semi-degenerate core a few years prior to a SN explosion, and the 
strongest flash could eject most of the H envelope with velocities of 
~ 100 kms -1 . The wind speed inferred from the narrow compo- 
nent of Ha (about 200 kms -1 ), which has likely CSM as origin, 
has a larger velocity width than is seen normally in AGB stars, 
although such velocities are not unprecedented in post-AGB stars. 
Super-AGB stars, or ECSN progenitors, may be distinguished from 
very massive progenitors of similar luminosity by their much cooler 
effective temperatu res ( < 3000 K for ECSN vs. ~ 3400 K for 
CCSN progenitors dWanaio et al] |2008)). The exciting star in the 
DUSTY model described above has a t emperature of 3000 K , hence 
consistent with cool super-AGB stars (Eldridge et al]2007T) . A lack 
of q elements in the case of E CSN is also a key to distinguish 
betwe en the ECSN and CCSNe dKitaura et afll2006tlWanaio et alj 
2008). It is interesting to note that the Crab remnant has also been 
suggested to have originated in an ECSN g iven the low kinetic 
energy and a sm all amount of a elemen ts dNomotoetal.il 19821 ; 
IChevaliei| [l984). Moreover ISwartz et all dl99lh . as already dis- 
cussed in the previous section, suggested that the ECSN progenitor 
may lose a large fraction of their envelope and become a type II-L 
SNe. 

Finally the rate of these types of transients would be a use- 
ful guide to their physical origins. IWanaio et alj (2008) suggest, 



from nucleosynthesis argu ments, that ECSNe m ust be ^ 30% of 
all CCSN events, while iPoelarends et alj ( I2008h suggest that they 
should be s omewhat rarer (~ 7-8 %) if they come from a narrow 
mass range. iThompson et alj d2008h propose that there are possibly 
4 known events that are similar, and that their faint peak magnitudes 
mean that many more may be undis covered in th e Loca l Universe. 
The local SN sample compiled by ISmartt et alj (2009) shows 92 
CCSNe within a local 28Mpc volume, within a 10.5 yr period. At 
le ast there are 4 candidate s which are SN2008S-like, as discussed 
in lThompson et al] d2008h . Hence a lower limit to the rate of these 
events of ^ 5% can be inferred, which is not inconsistent with the 
arguments from theory. 



8 CONCLUSIONS 

In summary, the combination of our monitoring data and the evi- 
dence from the progenitor studies suggests that a weak SN explo- 
sion, due to core-collapse through the electron-capture SN mech- 
anism, of a massive star with an initial mass around 6-10 Mq 
in the super-AGB stage is a plausible explanation for SN 2008S. 
By implication NGC 300 OT2008-1 and M 85 OT2006-1 are 
likely of similar origin. The progenitor star was not a normal RSG 
which should produce the standard type II-P SNe ( Smartt et al. 
2009). However, if the pre-explosion MIR flux reflects the stel- 
lar photospheric flux, the object has a luminosity similar to mod- 
els of super-AGB stars. Extremely high mass-loss rates could cre- 
ate cir cumstellar shells whi ch are optically thick, cool, dense and 
dusty ( IThompson et al.l l2008h . Such high mass-loss through ther- 
mal pulses has been predicted by theory and ob scured massive 
oxyge n-rich AGB stars are certainly observed by Ivan Loon et al.1 
J2005h 

The evidence that SN 2008S was indeed a SN explosion 
comes from two observations. Firstly, SN 2008S is similar in the 
total radiated energy to other faint SNe, and it shows moderate 
velocities of about 3000 kms -1 . The latter are hard to reconcile 
with either an LBV-like or non-destructive stellar eruption. Sec- 
ondly, we detect a tail phase which has a decay slope matching 
that of radioactive 56 Co. We estimate the mass of 56 Ni ejected to 
be 0.0014 ± 0.0003 M Q , which is marginally lower than that ob- 
served in the faintes t known type II - P SNe (^0.002-0.008 Mflfrom 
iTuratto et al] dl998l) ; lBenetti et alj d20oTh ; |Pastorello et alj d2004l) ) 
but is close enough that SN 2008S could be an extension of these 
low- energy explos i ons. 

IWanaio et al] d2008l) have presented models of electron cap- 
ture supernovae of progenitor AGB stars with an O-Ne-Mg core 
and an initial mass of 8.8Mgthat are in plausible agreement with 
SN 2008S an d similar tran s ients. These weak explosions, as also 
discussed by IKitaura et al] d2006l) . give an ejected 56 Ni mass of 
0.002-0.004 M Q . SN 2008S is also significant in that, apart from 
SN 1987A, it has allowed the earliest ever MIR observation of a 
SN, and the remarkably large flux detected confirms the presence 
of substantial circumstellar material around the progenitor star. In 
addition, SN 2008S developed a NIR excess at later times. This can 
be interpreted as optically thin thermal emission from ~ 10 -6 M 
of amorphous carbon grains or ~ 10 - Mq of silicate grains. How- 
ever, larger masses of grains, which would be optically thick, are 
not ruled out. The grain location and heating mechanisms are un- 
certain. While some of the NIR emission may originate in new, 
radioactively-heated ejecta grains, the magnitude of the later NIR 
flux requires that at least a proportion must come from an additional 
source. This might be reverse shock heating of new ejecta grains. 
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Alternatively it might be due to shock heating of new (CDS) dust 
or old dust in the circumstellar region. 

Our suggested scenario can be tested in the future. If the star 
has cataclysmically exploded as a SN, then, once the remnant fades, 
there should be no source left with a luminosity similar to the pro- 
genitor (10 4,6 Lq). This needs to be tested over the entire optical 
to MIR wavelength range to ensure that no progenitor star remains, 
whether exposed or concealed by newly-formed dust. In addition, 
if the SN ejecta become visible, one might expect to see broad lines 
of forbidden oxygen (e.g. [O I] AA 6300,6364) and other interme- 
diate mass elements. This may be difficult, as the events are in- 
trinsically faint, but it could be possible in particular for NGC 300 
OT2008- 1 as its distance is only 2 Mpc and hence it can be mon- 
itored for a long time. The tail phase, following the 56 Co decay, 
should continue at least for another ~ 300 days, and the slope of 
this decay should be monitored closely in both SN 2008S and NGC 
300 OT2008-1. Finally the true rate of these transients will give 
further insights into their nature, and deeper searches by future sky 
surveys may discover more of these events dYoung et alj2008h . 
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Table Al. Magnitudes of the local sequence stars in the field of SN2008S. 



ID 


U 


B 


V 


R 


I 


1 


16.41 ± 0.04 


16.29 ± 0.04 


15.45 ± 0.02 


14.97 ± 0.02 


14.46 ± 0.05 


2 


14.56 ± 0.02 


14.29 ± 0.02 


13.55 ± 0.02 


13.13 ± 0.01 


12.77 ± 0.02 


3 


18.35 ± 0.04 


17.85 ± 0.03 


16.85 ± 0.02 


16.30 ± 0.03 


15.75 ± 0.03 


4 


14.77 ± 0.02 


14.50 ± 0.02 


13.76 ± 0.01 


13.34 ± 0.02 


12.95 ± 0.03 


6 


16.83 ± 0.03 


16.61 ± 0.03 


15.80 ± 0.02 


15.32 ± 0.03 


14.83 ± 0.02 


7 


17.94 ± 0.04 


17.62 ± 0.04 


16.68 ± 0.03 


16.16 ± 0.02 


15.58 ± 0.02 


12 


16.79 ± 0.03 


16.60 ± 0.02 


15.83 ± 0.02 


15.37 ± 0.02 


14.89 ± 0.01 


13 


18.41 ± 0.05 


17.57 ± 0.04 


16.43 ± 0.03 


15.75 ± 0.02 


15.14 ± 0.02 


16 


15.87 ± 0.02 


15.59 ± 0.02 


14.74 ± 0.02 


14.27 ± 0.03 


13.81 ±0.03 


17 


17.88 ± 0.03 


17.71 ±0.04 


16.91 ± 0.04 


16.45 ± 0.03 


15.97 ± 0.03 


19 


16.13 ± 0.02 


15.78 ± 0.02 


14.85 ± 0.02 


14.34 ± 0.02 


13.88 ±0.02 


20 


19.89 ± 0.06 


18.12 ±0.05 


16.86 ± 0.03 


15.99 ± 0.03 


15.32 ± 0.03 


21 


16.75 ± 0.04 


16.37 ± 0.03 


15.49 ± 0.02 


14.99 ± 0.02 


14.53 ± 0.02 


22 


16.55 ± 0.03 


16.44 ± 0.03 


15.66 ± 0.02 


15.21 ± 0.02 


14.74 ± 0.02 


23 




15.56 ± 0.02 


13.82 ± 0.03 


13.35 ± 0.03 


12.91 ± 0.03 


26 




14.55 ± 0.03 


14.38 ± 0.02 


13.70 ± 0.02 


13.09 ±0.02 
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Table A2. Emission features observed on SN 2008S spectra. 



Identification Rest Wavelength Observed Wavelength Notes 

(A) (A) 





4101 


4103 


blend 


Fell (28) 


4178.85 


4182 


blend 


Fell (27) 


4385.38 


4387 




Fell (37) 


4472.92 


4475 




Fell (37) 


4491 


4493 




Fell (38) 


4508.28 


4511 


blend 


Fell (37) 


4520 


4522 




Fell (38) 


4541.52 


4543 


blend 


Fell (38) 


4549.47 


4553 




Fell (38) 


4576.331 


4578 




Fell (38) 


4583.829 


4586 




Fell (38) 


4620.5 


4622 




Fell (186) 


4635.3 


4637 




Bp 


4861 


4867 


blend 


Fell (42) 


4923.92 


4926 




Fell (42) 


5018.43 


5020 




Fell (42) 


5169.0 


5172 




Fell (49) 


5197.56 


5199 




Fell (49) 


5234 


5237 




Fell (49) 


5254.9 


5257 




Fell (48) 


5264.8 


5267 




Fell (49) 


5275.99 


5277 




Fell (49) 


5316.6 


5318 




Fell (49) 


5325.5 


5328 




Fell (48) 


5337.71 


5340 




Fell (49) 


5425.3 


5428 




Nal 


5890 5896 


5892 




Fell (46) 


5991.38 


5995 




Fell (46) 


6084.11 


6089 




Fell (46) 


6113.33 


6116 




Fell (74) 


6238.38-6239.9 


6243 


blend 


Fell (74) + Fell (92) 


6247.55 -6248.89 


6250 


blend 


[OI] 


6300.23 


6304 




[OI] 


6363.88 


6368 




Fell (40)+ Fell (93) 


6369.5-6371 


6372 




Fell (74) 


6416.9 


6418 




Fell (40)+ Fell (199) 


6432.7- 6433.9 


6433-6440 


blend 


Fell (74) 


6456.38 


6456-6461 


blend 


Fell (92)+ Fell (92) 


6491.2 + 6493.035 


6492-6496 


blend 


Fell (40)+ Fell (92) 


6516.1 +6517.0 


6519 


blend 


H a 


6563.5 


6568 




[FeII](U) 


7155.14 


7159 


? 


Fell (73) 


7222.39 


7226 




Fell (73) 


7224.51 


7228 




[Call] 


7291 


7295 




Fell (73) 


7310.24 


7313 




[Call] 


7324 


7327 




Mnll (4) 


7353.52 


7358 


? 


Mnll (4) 


7415.80 


7420 


? 


Mnll (4) 


/43Z.Z / 


n A 1 Q 

/43o 




Fell (73) 


7449.3 


7453 


blend 


Fell (73) 


7462.38 


7465 




Fell (73) 


7515.9 


7520 




Fell (72) 


7533.42 


7538 




Fell (73) 


7711.71 


7716 




Call (2) 


8498 


8505 




Call (2) 


8542 


8547 




Call (2) 


8662.14 


8669 




01(4) 


8446.35 




ph 256 
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